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PREFACE 


This  report  documents  the  work  performed  by  Computer  Sciences  Corporation  under 
Contract  DAAH03-75-A-0045.  The  task  reported  here  is  concerned  with  signal 
processing  for  a  phase  coded  CW  radar  being  developed  at  the  Advanced  Sensors 
Directorate  of  the  U.S.  Army  Missile  Research  and  Development  Command, 

Redstone  Arsenal,  Alabama.  Computer  Sciences  Corporation  wishes  to  take  this 
opportunity  to  express  its  appreciation  to  those  that  have  contributed  directly 
and  indirectly  to  the  work  reported  herein,  in  particular,  Mr.  W.  L.  Low, 

Dr.  D.  W.  Burlage  and  Mr.  R.  R.  Boothe  of  the  Advanced  Sensors  Directorate. 
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SECTION  1  -  INTRODUCTION 


The  need  for  a  short  range,  low  altitude  air  defense  system  that  is  capable  of 
operating  at  all  times  and  in  adverse  weather  conditions  dictates  the  use  of 
RF  devices.  Passive  RF  systems  which  depend  entirely  on  emissions  from 
attacking  aircraft  can  easily  be  rendered  useless  through  emission  turnoff  by 
the  threat  aircraft.  Consequently,  air  defense  systems  have  the  need  for  an 
active  RF  capability,  such  as  a  radar,  in  order  to  fulfill  their  task. 

However,  with  the  advent  of  antiradiation  missiles  (ARM)  the  susceptibility 
of  active  radars  to  detection  by  electronic  intelligence  receivers  and  attack 
radar  warning  receivers  has  become  a  major  problem.  Survivability  of  such  a 
radar  can  be  increased  by  designing  it  so  as  to  deny  detection  and  classification 
by  the  ARM  receiver. 

The  detection  performance  of  a  radar  depends  on  the  signal  energy  or  the 
average  power.  Thus  for  a  given  detection  range,  and  hence  average  power,  the 
peak  power  can  be  reduced  by  increasing  the  pulse  width,  and  continuous  wave 
(CW)  transmission  requires  the  minimum  peak  power.  This  reduces  the  output 
power  requirements  of  the  transmitter;  and,  more  importantly,  the  radar 
detectability  by  broadband  ARM  receivers  is  reduced.  In  addition,  CW 
transmission  complicates  the  sorting  and  designation  process  by  the  ARM 
receiver  in  that  the  commonly  used  discriminants  of  pulse  repetition  frequency 
and  pulse  width  are  not  available.  One  other  significant  characteristic  of 
CW  transmission  is  the  inherent  multipath  problem  generated  by  such  signals. 

This  is  due  to  the  fact  that  CW  signals  have  no  leading  edges  while  ARM  receivers 
rely  heavily  on  leading  edge  gating  to  reduce  multipath  effects.  The  multipath 
problem  for  the  receiver  will  be  more  severe  if  the  radar  employs  techniques  such 
as  intentionally  illuminating  the  ground  or  other  multipath  scatterers. 


A  high  gain  pencil  beam  antenna  should  be  used  to  provide  a  sufficient  power 
density  at  the  target  with  a  low  transmitter  power  level.  This  low  power 
coupled  with  low  sidelobes  on  transmit  provides  an  extremely  low  detectability 
by  a  sidelobe  ARM.  Moreover,  since  the  radar  can  be  detected  in  the  main  beam 
during  surveillance,  a  small  beam  size  confines  this  detectability  to  a  small 
portion  of  the  total  coverage. 

When  a  battery  of  radars  is  operating  in  the  area,  frequency  agility  can  be 
used  to  complicate  the  designation  process  in  the  ARM  receiver.  Since  mainlobe 
detection  of  the  rader  is  possible,  beam-to-beam  frequency  agility  makes  this 
detection  essentially  useless  when  the  radar  is  not  an  isolated  radar.  The 
possibility  of  spot  jamming  by  standoff  jammers  is  virtually  eliminated  by 
frequency  agility. 

Since  a  CW  signal  without  modulation  does  not  provide  range  resolution,  it  is 
necessary  to  increase  the  bandwidth  of  the  signal  by  using  a  suitable 
modulation.  Selecting  the  type  of  waveform,  and  hence  the  bandwidth,  is  a 
compromise  between  potentially  conflicting  performance  criteria  such  as  high 
resolution,  unambiguous  range  and  doppler  measurements,  adequate  performance  in 
clutter  and  ECM  environments,  low  power  density  for  increased  covertness, 
minimum  complexity  and  cost  of  waveform  generation,  signal  processing,  etc. 
Sawtooth  frequency  modulation  and  a  binary  phase  modulated  CW  carrier  are 
prime  candidates  for  the  radar  waveform. 

The  fundamental  disadvantages  of  a  phase  coded  CW  radar  are  inherent  to  CW 
radars:  antenna  isolation  and  large  dynamic  range  caused  by  close-in  clutter 
and  antenna  spillover.  This  report  documents  a  study  of  signal  processors  for 
a  binary  phase  coded  CW  radar,  where  a  maximum  length  shift  register  generated 
sequence  is  used  for  phase  modulation.  Following  a  brief  review  of  phase 
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coded  CW  radars  and  the  problems  associated  with  such  radars,  several  analog, 
digital  and  hybrid  signal  processor  configurations  are  proposed.  The 
advantages  and  drawbacks  of  each  configuration  are  discussed.  In  order  to 
evaluate  the  performance  of  the  proposed  processors,  digital  computer 
simulations  of  the  input  signal  and  the  processors  are  developed.  These 
simulations  are  used  to  demonstrate  the  feasibility  of  digital  processing  in 
spite  of  the  large  dynamic  range. 
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SECTION  2  -  PHASE  CODED  CW  RADAR 


2.1  INTRODUCTION 

Since  an  unmodulated  CW  radar  cannot  resolve  targets  in  range,  some  form  of 
modulation  is  required.  The  ability  of  the  radar  system  to  resolve  targets  in 
range  depends  on  the  bandwidth  of  the  transmitted  modulation;  in  fact  the  range 
resolution  is  inversely  proportional  to  the  bandwidth.  Expansion  of  bandwidth 
is  achieved  by  modulating  the  CW  carrier  with  a  repetitive  modulation  waveform. 
Two  such  waveforms  which  are  most  commonly  used  are  the  sawtooth  frequency 
modulation  and  binary  phase  coding.  In  addition  to  improved  range  resolution, 
the  enlarged  bandwidth  produces  a  lower  power  density  thereby  increasing 
covertness  against  ARM  systems  that  use  narrow-band  superhet  receivers. 


2.2  BINARY  PHASE  CODED  SIGNALS 

The  simplest  form  of  phase  coding  is  the  binary  or  the  phase-reversal  code  in 

o 

which  the  phase  of  the  carrier  is  shifted  by  either  zero  or  180  .  The  complex 
representation  of  such  a  signal  is  given  by, 

i/i(t)  =  u(t)  e^  °  (2.1) 


where  w0  is  the  carrier  frequency  in  radians/sec  and  u(t)  is  the  modulation 
waveform.  For  a  binary  phased  coded  signal,  the  modulation  waveform  is  of  the 
type 


-1* 

u(t)  =  e  n  for  (n-l)6  <t  n6 
$n  =  0  or  180° 

where  <5  is  the  code  clock  period. 


1,2,, 


(2.2) 
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From  (2.2)  it  is  clear  that  u(t)  is  a  rectangular  wave  assuming  values  of  +1  or 
-1,  and  can  be  written  as 


u<c>  “  XJ  an  Vc> 

where 

^(t)  =  1  for  (n-l)8  <t  <_  n<$ 

=  0  elsewhere 


(2.3) 


(2.4) 


and 

an  =  ±1  .  (2.5) 

An  example  of  a  binary  modulating  waveform  and  the  corresponding  modulated  carrier 
are  shown  in  Figure  2.1.  Selection  of  the  modulating  waveform  is  equivalent  to 
choosing  the  coefficients  a^,  ....  so  that  the  signal  has  the  desired  range 

resolution  properties. 

2.3  PSEUDO-NOISE  SEQUENCES 


Among  the  most  commonly  used  binary  modulating  sequences  are  the  Barker  codes  and 
the  pseudo-noise  (PN)  sequences,  both  of  which  are  periodic  and  have  similar  range 
resolution  properties.  Barker  codes  suffer  from  a  major  disadvantage  in  that 
their  existence  has  been  established  only  for  codes  with  short  periods.  Pseudo¬ 
noise  sequences,  on  the  other  hand,  can  be  easily  generated  using  linear  shift 
register  generators  and  have  been  proved  to  exist  for  all  lengths  of  the  type 
N  =  2n-l.  Although  pseudo-noise  sequences  exist  for  numerous  other  lengths,  only 
the  linear  maximum  length  shift  register  generated  sequences  or  so-called  linear 
m-sequences  are  considered  in  this  study. 
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Linear  maximum  length  sequences  can  be  obtained  from  the  system  shown  In 
Figure  2.2.  This  shift  register  generator  of  degree  n  consists  of  n  storage 
elements  (n  stages),  a  clock  pulse  generator  which  generates  pulses  at 
intervals  of  s,  and  a  feedback  logic  circuit.  When  a  clock  pulse  arrives,  the 
contents  of  each  stage  are  transferred  to  the  next  stage,  and  the  content  of 
the  last  stage  is  the  output.  The  new  state  of  the  first  stage  is  a  linear 
Boolean  function  of  the  previous  contents  of  some  or  all  the  stages.  For  each 
stage,  the  OFF  state  is  symbolized  by  0  and  the  ON  state  by  1.  The  binary 
output  of  the  feedback  logic  unit  is  of  the  form 

f(xj_,  x2,  .  .  .  ,  xn)  ■=  ©C2x2@.  .  .©C^  (2.6) 

where  x^  denotes  the  state  of  the  i1"*1  stage,  each  constant  is  either  zero  or 
one  and  the  symbol  ©  denotes  modulo  2  addition.  It  can  be  easily  shown  that 
the  output  sequence  of  this  generator  is  periodic  whose  period  is  at  most 
N  =  2n-l.  Of  the  total  of  2n  possible  feedback  combinations,  only  a  few  result 
in  sequences  whose  period  is  N  and  such  sequences  are  called  linear  maximal 
length  sequences.  The  feedback  connections  which  yield  maximal  length  sequences 
have  been  extensively  tabulated  and  a  part  of  it  is  shown  in  Table  2.1.  The  use 
of  this  table  will  be  demonstrated  with  an  example.  Consider  a  6-stage  generator, 
n  =  6  for  which  the  maximum  period  is  N  =  63.  The  table  lists  three  possible 
feedback  connections: 

(i)  [6,1]  which  corresponds  to 

f(xi,  x2,  x3,  x4,  x5,  x6)  =  x6©  xi  (2.7) 

(ii)  [6,5, 2,1]  which  corresponds  to 

f(xx,  x2,  X3,  x4,  X5,  X6)  =  X6©  X5  ©  X2  ©x3  (2.8) 
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Table  2.1  Feedback  Connections  for  Linear  m-sequences 


Number  of 
stages 


Code 

length 


Feedback  Connections 


12,1] 

13.1) 

14.1) 

[5,2]  15,4,3,21  [5, A, 2,1] 

(6.1]  [6,5,2,1]  [6,5,3,2] 

[7,11  17,3]  [7,3,2,1]  [7, A, 3, 2] 

[7 ,6, A ,2]  [7, 6, 3,1]  [7, 6, 5, 2] 

(7,6, 5, A, 2,1]  [7, 5, A, 3, 2,1] 

[3, A, 3, 2)  [8,6,5,31  [8, 6, 5, 2] 

(8.5.3.1]  [8,6,5,11  [8,7,6,11 

18.7.6.5.2.1)  [8, 6, A, 3, 2,1] 

[9, A)  [9, 6, A, 3)  [9, 8, 5, A]  [9,8,A,1] 

(9.5.3.2]  [9,8,6,51  [9, 8, 7, 2] 

[9, 6, 5, A, 2,1)  [9, 7, 6, A, 3,1] 
[9,8,7,6,5,31 

[10.3]  [10,8,3,2]  [10, A, 3,1]  [10,8,5,1] 
[10, 8, 5, A]  [10,9,A,1]  [10, 8, A, 3] 

(10.5.3.2]  [10,5,2,1]  [10, 9, A, 2] 

[11.1]  [11,8,5,2]  [11,7,3,2]  [11,5,3,2] 

(11.10.3.2]  [11,6,5,1]  [11,5,3,1] 

[11, 9, A, 11  [11,8,6,2]  [11,9,8,31 

[12,6,A,1]  [12,9,3,2]  [12,11,10,5,2,1] 
[12, 11, 6, A, 2,1]  [12,11,9,7,6,5] 

[12.11.9.5.3.1]  [12, 11, 9, 8, 7, AJ 
[12,11,9,7,6,5]  [12,9,8,3,2,1] 

[12.10.9.8.6.2] 

(13, A, 3,1]  [13, 10, 9, 7, 5, A] 

(13, 11, 8, 7, A, 1]  [13,12,8,7,6,5] 

[13.9.8.7.5.1]  [13, 12, 6, 5, A, 3) 

[13.12.11.9.5.3]  [13,12,11,5,2,1] 

(13, 12, 9, 8, A, 2]  [13, 8, 7, A, 3, 2) 

[14.12.2.1]  [1A.13.A.2]  [14,13,11.91 
(1A,10,6,1]  IlA.11,6,1]  [1A, 12, 11,1] 

[ 1  A, 6, A, 2]  [1A, 11, 9, 6, 5, 2) 


(iii)  [6,5, 3, 2]  which  corresponds  to 


x2,  x3,  x4,  x5,  xfi)  -  x6  ©  x5  ©  x3  ©  x2  .  (2.9) 

For  a  given  feedback  connection,  different  initial  conditions  in  the  shift  registers 
produce  the  same  sequence  with  different  shifts.  The  initial  condition  of  all  zeros 
produces  a  sequence  of  all  zeros  and  is  therefore  prohibited.  The  output  of  the 
shift  register  generator  is  a  sequence  of  zeros  and  ones  which  correspond  to  phase 
shifts  of  zero  (a^  =  +1)  and  180°  (a^  =  -1). 


The  fundamental  properties  of  a  linear  m-sequence  produced  by  an  n-stage  generator 
are: 

(i)  the  sequence  is  periodic  with  period  N  =  2n-l 

(ii)  in  each  period,  the  difference  between  the  number  of  ones  and  the 
number  of  minus  ones  is  one 


(iii)  if  an  m-sequence  is  multiplied  by  a  shifted  version  of  itself,  the  resulting 
sequence  is  also  a  shifted  version  of  the  original  sequence 


(iv)  the  periodic  autocorrelation  is  two-valued,  i.e., 
N 

ai  ai+k  *  n  if  1=0 *  ±n»  ±2n»  •  •  • 

Val 

=  -1  otherwise  . 


(2.10) 


2.4  PERFORMANCE  OF  PHASE  CODED  CW  RADARS 

The  ability  of  the  radar  system  to  resolve  targets  in  range  depends  on  the 
autocorrelation  of  the  transmitted  modulation.  The  ideal  autocorrelation  has  a 
peak  at  zero  delay  and  is  zero  elsewhere.  Consider  a  PN  sequence  of  period  N  and 
clock  period  6.  The  code  repetition  interval  is  A  ■  N6.  The  normalized 
autocorrelation  of  such  a  sequence  can  be  shown  to  be  periodic  with  period  A  and 


has  a  shape  as  shown  in  Figure  2.3.  Due  to  the  periodicity  of  the  autocorrelation, 
the  power  spectrum  is  a  line  spectrum  and  it  can  be  shown  that  its  envelope  is 


given  by 


U(f) 


5jiL 3 
7L 


=  A  sa2  (irf/B)  (2.11) 

where  B  **  l/<5,  except  at  DC  where  the  value  is  A/N^.  The  magnitude  spectrum  of  the 
waveform  (square  root  of  the  power  spectrum)  is  sketched  in  Figure  2.4.  From  the 
magnitude  spectrum,  it  can  be  seen  that  the  bandwidth  of  the  signal  is  approximately 
equal  to  B,  the  basic  clock  rate. 


While  the  use  of  a  veiy  wide  bandwidth  signal  offers  some  improved  covertness 
against  ARM  receivers  which  use  narrowband  superhet  receivers,  it  also  has  several 
disadvantages  in  regard  to  the  radar  performance.  The  very  high  range  resolution 
caused  by  the  increased  bandwidth  poses  such  problems  as  range-gate  fly through 
and  the  need  for  more  receiver  channels  to  cover  all  the  desired  range  cells  in  the 
required  time.  A  wide  bandwidth  signal  is  more  likely  to  encounter  spot  jammers 
spaced  randomly  across  the  operational  band.  The  number  of  frequency  bands 
available  for  frequency  agility  is  sharply  reduced  by  the  use  of  wide  bandwidth 
signals.  Finally,  the  large  bandwidth  will  have  an  undesirable  cost  impact  on  the 
RF  portion  of  the  radar  and,  in  particular,  the  low  sidelobe  antenna.  Although 
very  large  bandwidths  are  not  desirable,  the  bandwidth  must  be  sufficiently  large 
so  as  to  provide  good  range  resolution  and  improved  performance  in  heavy  clutter. 

The  substantial  benefits  of  a  modulated  CW  radar  are  somewhat  tempered  by  the 
inherent  drawbacks  of  a  CW  radar.  Chief  among  these  are  the  multipath  effects, 
antenna  spillover  and  large  receiver  dynamic  range  requirements  caused  by  close-in 
clutter.  Antenna  spillover  can  be  minimized  by  properly  shielding  the  receive 
antenna  from  the  transmit  antenna.  The  task  of  designing  a  signal  processor, 
especially  a  digital  processor,  is  made  considerably  more  difficult  by  the  large 


dynamic  range  requirement. 
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2.5  TYPICAL  PARAMETER  VALUES 


The  following  is  a  list  of  symbols  and  their  typical  values  that  will  be  used 
throughout  the  report  to  demonstrate  the  operation  of  the  signal  processors. 


Carrier  Frequency: 

fo 

= 

10  GHz 

Wavelength: 

X 

= 

0.03  Meters 

Intermediate  Frequency: 

f . 

i 

= 

30  MHz 

Clock  Rate: 

B 

= 

5  MHz 

Range  Resolution: 

AR 

= 

30  Meters 

Code  Length: 

N 

= 

63 

No.  of  Shift  Register  Stages: 

n 

= 

6 

Unambiguous  Range: 

®unamb 

=  1890  Meters 

Code  Repetition  Period: 

A 

s= 

12.6  n  sec 

Code  Repetition  Frequency: 

*r 

£ 

80  KHz 

Look  Time: 

T 

= 

2  m  sec 

No.  of  Code  Periods  in  One  Look: 

K 

= 

158 

Maximum  Range  Rate  of  Interest: 

v 

vmax 

=  360  meters/sec 

Maximum  Doppler: 

fd 

=  24  KHz 

Note  that  due  to  the  relatively  small  unambiguous  range,  additional  means  such  as 
varying  the  pulse  repetition  rate  must  be  employed  to  resolve  range  ambiguities. 

The  existence  of  the  dynamic  range  problem  can  be  easily  demonstrated  by  considering 
a  1  target  at  15  KM  range  and  a  10^  m^  fixed  clutter  at  a  2  KM  range.  For  such 
a  situation  the  signal-to-clutter  ratio  is 


-75  dB 
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In  addition  to  antenna  spillover  and  noise,  various  types  of  clutter  must  be  con¬ 
sidered  in  the  design  of  a  signal  processor.  These  include  in-range  clutter  (in  the 
same  range  cell  as  the  target),  out-of-range  distributed  clutter,  fixed  clutter  such 
as  buildings  and  weather  clutter. 


SECTION  3  -  THE  ANALOG  PROCESSOR 


3.1  INTRODUCTION 

The  signal  processor  described  in  this  chapter  is  a  simple  analog  processor  and 
is  presented  to  demonstrate  the  principles  of  signal  processing  for  phase  coded 
CW  radars.  Some  of  the  problems  associated  with  such  a  system  are  also  elaborated. 
Typical  parameter  values  set  forth  in  the  previous  chapter  are  used  to  explain  the 
operation  of  the  processor.  It  must  be  noted  that  this  processor  is  similar  to  the 
one  being  used  to  demc.istrate  the  feasibility  of  the  phase  coded  CW  radar  concept. 

3.2  DESCRIPTION  OF  THE  PROCESSOR 

A  simplified  schematic  of  the  processor  is  shown  in  Figure  3.1.  The  received  RF 
signal  is  mixed  down  to  a  convenient  IF  frequency  f^  and  amplified.  This  signal  is 
then  sent  into  63  parallel  range  channels  through  a  power  divider.  Each  channel 
contains  a  decoder  or  a  code  demodulator  where  the  incoming  signal  is  multiplied 
by  the  binary  code  but  with  a  different  shift  in  each  channel.  The  target  return 
may  be  written  in  the  form, 

gT(t)  =  p(t-t)  cos  [2xt(fi+fd)  +  <f>]  (3.1) 

where  p(t)  is  the  periodic  code, 

fd  is  the  target  doppler  and 
t  is  the  range  delay. 

The  in-range  clutter  has  a  similar  form, 

gc  (t)  “  a(t)  p(t-T)  cos  (2irtfi  +  'J'c)  (3.2) 

where  a(t)  is  the  low  frequency  amplitude  variation. 


Figure  3.1  The  Analog  Processor 
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The  out-of-range  distributed  clutter  is  the  sum  of  the  returns  from  several  range 
cells,  each  of  the  form 

gd  (t)  =  b  (t)  p(t-Td)  cos  (2tftfi  +  <j>d)  (3.3) 

where  t  ,  t  x . 
d 

One  of  the  range  channels  contains  a  decoder  which  has  the  same  delay  as  the 
target.  This  channel  will  be  called  the  in-range  channel  and  the  other  channels 
the  out-cf-range  channels.  The  decoder  output  due  to  the  target  in  the  in-range 
channel  is, 

rT(t)  =  gT(t)  p(t-x) 

=  p(t-x)  p(t-x)  COS  [2rrt  (f-i+fd)  +  ] 

=  cos  [2irt(f^+fd)  +  <f>]  (3.4) 

since  p(t-x)  can  only  assume  values  +1  and  -1. 

The  decoder  output  in  an  out-of-range  channel  due  to  the  target  is, 

ST(t)  =  gT(t)  p(t-T1)  ,  7*  T 

=  p ( t— x)  p(t-X2)  COS  [27Tt(fi+fd)  +  |J> ) 

=  p(t-x2)  cos  [2nt(fi+fd)  +  <|>]  (3.5) 

where  x2  ^  0  and  use  is  made  of  the  fact  that  a  PN  sequence  multiplied  by  its 
shifted  version  yields  the  same  sequence  with  a  different  shift.  The  output  of 
the  decoder  is  filtered  by  a  bandpass  filter  with  a  clutter  notch  centered  at  f^. 

Let  the  outputs  of  the  bandpass  filter  be  uj(t)  and  vT(t)  when  the  inputs  are  rT(t) 
and  Sj(t)  respectively.  The  spectra  of  the  signals  in  the  in-range  and  out-of-range 
channels  due  to  the  target  return  are  shown  in  Figures  3.2  and  3.3. 


3-3 


(c)  Filter  Characteristics 


Figure  3.2  Target  Return:  In-Range  Channel 


(b)  Spectrum  After  Decoding 


(d)  Output  Signal  Spectrum 


Figure  3.3  Target  Return:  Out-of-Range  Channel 


Note  that  the  line  spectrum  of  Figure  2.4  has  now  been  replaced  by  a  spectrum 
of  non-zero  width.  This  is  due  to  the  finite  look  time  of  2  msec  which 
corresponds  to  a  spectral  width  of  500  H2. 

The  decoder  outputs  due  to  the  in-range  and  out-of-range  clutter  can  be  written 
as, 

rc  (t)  =  a(t)  cos  (2ntfi  +  (}>c)  (3.6) 

rd  (t)  *•  b (t)  p(t-T3)  cos  (2rrtfi  +  <j>d)  (3.7) 

where  T3  f  0.  Let  the  corresponding  outputs  of  the  bandpass  filter  be  uc  (t) 
and  uj  (t) .  The  spectra  of  the  clutter  signals  in  the  in-range  and  out-of-range 
channels  are  shown  in  Figures  3.4  and  3.5. 

The  output  of  the  bandpass  filter  is  passed  through  a  detector  and  then  integrated. 
The  output  of  the  integrator  is  presented  to  the  thresholding  device  for  detection. 

3.3  ANALYSIS  OF  THE  PROCESSOR 

Several  observations  can  be  made  from  an  examination  of  the  developments  of  the 
previous  section.  The  response  due  to  a  target  in  an  out-of-range  channel  is  nearly 
36  dB  below  that  in  the  in-range  channel.  Therefore,  the  effect  of  range  sidelobes 
is  negligible.  Improvement  in  signal-to-noise  ratio  is  derived  by  reducing  the 
bandwidth  from  the  original  5  MHz  to  50  KHz  at  the  bandpass  filter  output.  This 
represents  a  20  dB  gain  in  the  signal-to-noise  ratio.  Both  in-range  and  out-of¬ 
range  clutter  can  be  attenuated  by  the  notch  in  the  bandpass  filter.  The  notch 
is  6  KHz  wide  corresponding  to  a  doppler  frequency  range  of  -45  to  +45  meters/sec. 
Targets  with  radial  velocities  of  less  than  45  meters/sec  cannot  be  detected. 


gTjWfr 


36  dB 


(a)  Input  Clutter  Spectrum 


4o  dB 


(b)  Spectrum  After  Decoding 


3  KHz-H 


25  KHz^-j 


(c)  Filter  Characteristics 


Figure  3.4  In-Range  Clutter 
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(a)  Input  Clutter  Spectrum 


i 

l  3  KHzH^ 

l+25  KHz-H^ 

fi 

(c)  Filter  Characteristics 


Figure  3.5  Out-of-Range  Clutter 


The  dual  requirements  of  adequate  clutter  suppression  and  low  target  return 
attenuation  dictate  the  need  for  steep  transitions  at  the  notch.  This,  however, 
results  in  a  long  settling  time  for  inputs  in  the  stop  band.  Thus,  the  transient 
response  due  to  the  large  clutter  inputs  persist  for  a  long  period  of  time, 
thereby  severely  degrading  the  clutter  cancellation  performance  of  the  processor. 
The  problem  of  long  settling  time  may  be  alleviated  to  some  extent  by  using  a 
suitable  weighting  scheme  prior  to  filtering.  Even  with  the  weighting,  it  will 
be  necessary  to  wait  for  the  output  to  settle  before  starting  integration.  The 
weighting  and  the  loss  of  integration  time  result  in  a  reduced  gain  in 
signal-to-noise  ratio.  One  other  problem  associated  with  such  a  processor  is 
the  reinitialization  of  the  filter.  This  is  because  the  filter  must  be  cleared 
before  processing  begins  in  the  next  look  period.  The  problem  of  clearing  is 
not  straightforward  in  the  crystal  filters  needed  for  the  IF  bandpass  filtering. 

3.4  CONCLUSIONS 

A  simple  analog  signal  processor  was  presented  which  uses  several  range  channels, 
IF  notch  filtering  for  clutter  reduction,  and  non-coherent  integration.  The 
principle  of  operation  of  the  processor  was  explained  for  the  case  of  a  target  in 
clutter.  Eventhough  the  processor  is  effective,  it  suffers  from  the  disadvantage 
of  slowly  decaying  transient  response  and  the  problem  of  clearing  the  filter  prior 
to  the  next  look  period. 


SECTION  4  -  DIGITAL  PROCESSING 


4.1  INTRODUCTION 

As  explained  In  the  last  chapter,  one  of  the  major  drawbacks  of  the  analog 
processor  is  the  long  settling  time  in  the  clutter  rejection  filter.  An 
alternative  would  be  to  employ  digital  processing  where  it  is  much  easier  to 
design  filters  with  the  desired  characteristics.  For  example,  finite  impulse 
response  or  transversal  filters  can  be  designed  with  very  short  settling  times. 
Simple  examples  of  such  filters  are  the  commonly  used  two  and  three  pulse 
cancellers.  Other  advantages  of  digital  processing  over  analog  processing 
include  the  wide  flexibility  offered  by  digital  processing  and  the  possibility 
of  employing  coherent  integration  through  the  use  of  fast  Fourier  transform 
(FFT)  processors.  However,  digital  processing  has  one  significant  disadvantage 
compared  to  analog  processing,  it  being  the  additional  errors  introduced  due  to 
quantization.  This  chapter  discusses  some  of  the  factors  that  must  be  considered 
during  analog- to-digital  conversion  or  digitization. 

4.2  ANALOG-TO-DIGITAL  CONVERSION 

Analog- to-digital  (A/D)  conversion  consists  of  two  steps:  sampling  and 
quantization.  Sampling  is  the  process  by  which  the  analog  signal  is  converted 
into  a  discrete  time  signal.  The  time  interval  between  two  successive  samples, 
usually  constant,  is  called  the  sampling  interval  and  its  reciprocal  the  sampling 

rate  or  sampling  frequency.  While  the  sampled  signal  exists  only  at  discrete  instants 
of  time,  the  amplitude  can  assume  a  continuous  range  of  values.  Quantization 
is  the  process  where  the  amplitudes  are  approximated  by  a  discrete  set  of  values. 

This  is  necessary  because  in  digital  representation,  the  amplitude  must  be 
represented  by  a  finite  number  of  binary  bits  (finite  word  length).  For  example 


1024  discrete  levels 


if  the  word  length  is  10  bits,  then  there  are  only  2^  = 
of  amplitude  that  can  be  represented.  The  approximation  of  the  amplitude  by  a 
discrete  set  of  values  introduces  an  error  known  as  the  quantization  error 
or  quantization  noise. 

In  accordance  with  the  low-pass  sampling  theorem,  the  sampling  rate  must  be  at 
least  twice  the  highest  frequency  at  which  the  signal  has  a  significant  component. 
It  therefore  becomes  obvious  that  A/D  conversion  cannot  be  achieved  at  RF  or  at 
IF  because  of  the  unrealistic  sampling  rate  requirements.  The  signal  must 
therefore  be  shifted  down  to  video  prior  to  A/D  conversion.  Since  the  signals 
of  interest  are  narrowband  signals,  the  signal  can  be  brought  down  to  baseband, 
without  loss  of  information,  by  using  two  video  channels.  Hie  in-pliase  and 
quadrature  signals  can  be  obtained  from  the  IF  signal  by  simply  mixing  it  with 
sine  waves  at  IF  but  with  a  90°  phase  shift,  as  shown  in  Figure  4.1.  The 
in-phase  and  quadarture  signals  can  be  considered  to  be  the  real  and  imaginary 
parts  of  a  complex  video  signal.  The  spectrum  of  this  signal  is  similar  to  the 
one  shown  in  Figure  2.4  except  for  a  shift  due  to  the  doppler  frequency.  For  the 
set  of  typical  parameters  being  considered  in  this  study,  the  bandwidth  of  this 
signal  is  about  5  MHz.  Thus,  the  complex  video  signal  can  be  sampled  at  5  MHz 
without  appreciable  loss  oi  nformation. 

r  e  process  of  quantization  can  be  explained  with  a  simple  example.  Consider  an 
input  signal  whose  range  is  -1  to  1  volt  and  a  word  length  of  3  bits.  Assuming 
linear  quantization,  the  input-output  relationship  of  the  digitizer  is  as  shown  in 
Figure  4.2,  where  the  quantity  q  is  called  the  quantization  interval.  It  is 
evident  from  the  figure  that  if  z  is  the  middle  of  one  of  the  steps,  then  any 
input  in  the  range  z-q/2  to  z+q/2  yields  the  same  output.  The. effect  of 
quantization  can  therefore  be  looked  upon  as  the  introduction  of  an  error  whose 
value  varies  between  -q/2  and  q/2.  To  facilitate  the  analysis,  the  error  is 


Ill 


considered  to  be  a  random  variable  which  is  uniformly  distributed  between  -q/2 
and  q/2.  It  can  be  easily  shown  that  such  a  random  variable  has  zero  mean  and 
a  variance  given  by 


o2  =  q2/12  . 


(4.1) 


If  the  range  of  the  input  is  -V  to  +V  volts,  and  if  the  word  length  is  m-bits, 
then  it  can  be  easily  shown  that 


q  =  2V/2®  • 


(4.2) 


In  particular,  if  V  =  1  volt,  then 


_  2-(m-1) 


(4.3) 


The  dynamic  range  of  the  A/D  converter  is  defined  as  the  ratio  of  the  largest 

signal  power  to  the  quantization  error  power.  Assuming  a  sine  wave  input,  the 

largest  signal  that  can  be  converted  without  saturation  has  a  magnitude  of  V 
2 

and  a  power  of  V  /2.  Thus,  the  dynamic  range  in  dB  is, 


RdB  10  1Og10  q2/12 


6.02m  +  1. 76 


="  6m  dB 


(4.4) 


To  gain  insight  into  the  dynamic  range  problem  posed  by  the  CW  radar,  consider 
a  signal-to-clutter  ratio  of  -70  dB.  If  the  signal  power  is  1,  then  the  clutter 


variance  is 


and  clutter  standard  deviation  is 


(4.5) 


ac  =  3.16  x  10J  . 


(4.6) 


If  the  A/D  converter  input  range  is  adjusted  so  that  a  3ac  signal  is  not 
saturated,  then 

(4.7) 


(4.8) 

Thus,  the  total  excursion  of  the  target  return  is  less  than  q,  and  if  the  target 
return  were  not  corrupted  by  clutter  and  noise,  it  is  possible  that  the  output  of 
the  A/D  converter  is  a  constant,  and  that  the  target  return  will  be  completely 
lost.  One  way  to  avoid  this  problem  is  to  use  a  larger  word  length,  but  the  word 
length  is  limited  by  the  5  MHz  sampling  rate  requirement. 

4.3  PRELIMINARY  EXPERIMENTS 

Preliminary  simulation  experiments  were  performed  to  examine  the  feasibility  of 
detecting  a  sine  wave  in  clutter  when  the  signal  amplitude  is  considerably 
smaller  than  the  quantization  interval.  The  clutter  signal  included  both  DC 
clutter  and  random  AC  clutter  with  Gaussian  distribution  and  Gaussian  spectrum. 

A  sampling  rate  of  5  MHz  was  used,  and  a  matched  filter  matched  to  the  sine  wave 
for  a  period  of  2  msec  was  used  as  the  detector.  It  was  found  that  the  sine 
wave  could  be  detected  in  80  dB  clutter  with  a  word  length  as  small  as  10  bits. 
Based  on  these  results,  it  was  inferred  that  digital  processing  is  feasible  for 
the  CW  radar  despite  the  stringent  dynamic  range  requirements  posed  by  the  wide 
disparity  between  the  magnitudes  of  the  clutter  return  and  the  target  return. 


V  =  3a, 


=  9.5  x  10' 


With  a  13-bit  A/D  converter,  the  quantization  interval  is 


q  =  2V/213 

=  2.3  volts 
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4.4  OTHER  CONSIDERATIONS 

In  addition  to  the  quantization  errors  introduced  during  A/D  conversion,  roundoff 
errors  are  introduced  during  every  arithmetic  operation.  The  magnitudes  of  these 
errors  depend  on  the  word  length,  the  type  of  number  representation  (one's 
complement,  two's  complement,  etc.)  and  the  type  of  arithmetic  employed  (fixed 
point,  floating  point,  etc.).  In  the  case  of  a  complicated  processor,  such  as 
the  FFT  processor,  these  errors  are  substantial  and  must  be  taken  into  account  in 
the  design.  Another  effect  of  the  use  of  finite  word  length  is  the  change  in  the 
characteristics  of  a  filter  due  to  truncation  of  the  coefficients.  The  pole-zero 
structure  and  hence  the  frequency  response  are  altered  and  in  some  extreme  cases 
the  filter  may  become  unstable. 

It  is  not  necessary  to  maintain  a  constant  word  length  throughout  the  processor. 
The  word  length  and  the  quantization  interval  at  each  stage  must  be  chosen  so  as 
to  accomodate  the  ejected  dynamic  range  without  excessive  deterioration  due  to 
finite  word  length  effects.  The  word  length  at  each  location  in  the  processor 
must  be  minimized  from  the  standpoints  of  economy  and  speed  of  operation.  It 
is  therefore  of  fundamental  importance  to  design  the  processor  in  such  a  way 
that  the  dynamic  range  requirements  are  minimized. 
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SECTION  5  -  DIGITAL  PROCESSOR  CONFIGURATIONS 


5.1  INTRODUCTION 

Several  digital  processor  configurations  are  proposed  in  this  chapter.  The 
processors  are  based  on  two  different  concepts:  decoding  and  pulse  compression. 

The  analog  processor  described  in  the  previous  chapter  uses  the  decoding  approach. 

The  pulse  compression  method  is  based  on  converting  the  CW  signal  into  a  pulsed 
signal  followed  by  standard  pulse  radar  techniques.  Also  included  are  two  hybrid 
processors  where  some  analog  processing  is  performed  prior  to  A/D  conversion.  While 
reading  the  block  diagrams,  it  must  be  remembered  that  all  the  video  signals  are 
made  up  of  two  separate  signals,  that  is,  the  in-phase  and  the  quadrature  components 
of  the  original  high  frequency  narrowband  signal. 

5.2  PROCESSORS  USING  THE  DECODING  APPROACH 
5.2.1  Configuration  D-l 

A  schematic  of  this  processor  is  shown  in  Figure  5.1.  It  _an  be  easily  seen  that 
this  is  a  video  frequency  equivalent  of  the  analog  processor  discussed  in  Section  3, 
where  the  clutter  rejection  filter  has  been  replaced  by  a  low-pass  filter  with  a 
notch.  This  filter  serves  to  increase  the  signal-to-noise  ratio  since  it  passes  the 
desired  signal  while  reducing  the  noise  bandwidth  from  the  original  5  MHz  to  50  KHz. 
because  of  the  reduced  bandwidth,  the  output  of  the  filter  may  be  resampled  at  a  much 
lower  rate,  i.e.,  50  KHz,  without  loss  of  information.  Also,  since  clutter  is 
attenuated  by  the  notch  filter,  the  word  length  requirement  at  the  integrator  is  much 
lighter  than  that  at  the  A/D  converter. 

If  a  recursive  filter  is  used  for  the  notch/LP  filter,  it  suffers  from  the  same 
disadvantage  as  the  analog  filter  in  that  the  settling  time  is  very  large.  This 
t  n.  1 1  i  i  at  e  a  weighting  scheme,  before  the  filter  and  discarding  the  output 


Figure  5.1  Configuration  D-l 


until  the  transient  effects  are  negligible.  If  a  non-recursive  or  a  transversal 
filter  is  used,  the  transient  response  problem  is  almost  completely  eliminated. 

For  example,  if  a  L-tap  filter  is  used,  it  is  only  necessary  to  discard  the 
first  L  samples  of  the  output  to  avoid  transients.  This  is  usually  much  fewer 
than  the  number  of  samples  that  must  be  abandoned  when  a  recursive  filter  is 
used.  In  addition,  if  the  number  of  taps  is  the  same  as  the  ratio  of  the  input 
sampling  rate  to  the  reduced  sampling  rate  contemplated,  filtering  and  sampling 
rate  reduction  can  be  achieved  simultaneously  using  a  fixed  window  realization 
of  the  transveral  filter.  Such  a  filter  can  be  implemented  with  just  one 
multiplier  and  an  accumulator.  No  signal  weighting  is  necessary  when  non¬ 
recursive  filters  are  used.  Note  that  when  weighting  is  required,  it  can  be 
performed  prior  to  decoding.  Thus,  the  input  can  be  weighted  only  once  instead 
of  once  in  each  channel. 

5.2.2  Configuration  D-2 

As  can  be  seen  from  Figure  5.2,  this  is  the  same  as  D-l  except  that  non-coherent 
integration  is  replaced  by  coherent  integration.  This  is  achieved  by  using  an 
FFT  processor  with  128  input  samples,  and  a  doppler  resolution  of  500  Hz 
corresponding  to  2  msec  look  time.  The  obvious  advantage  of  coherent  processing 
is  the  additional  improvement  in  signal-to-noise  ratio  as  compared  to  the  non¬ 
coherent  case.  Another  advantage  is  the  doppler  resolution  made  available  by  the 
FFT  processor.  Even  if  doppler  measurement  is  not  contemplated,  doppler 
separation  may  be  useful  in  resolving  range  ambiguities.  It  is  shown  in 
Appendix  A  that  the  frequency  response  of  the  FFT  processors  feature  large  side 
lobes.  Thus,  if  there  are  two  targets  in  the  same  range  cell  but  with  different 
doppler  shifts,  the  presence  of  the  side  lobes  produces  a  large  interference  between 
the  two  target  returns.  Tt  Is  therefore  necessary  to  multiply  the  signal  by  a 
weighting  sequence,  before  fast  Fourier  transformation. 
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5.2.3  Configuration  D-3 

Since  the  FFT  processor  provides  frequency  resolution,  the  notch  filter  has  been 
removed  in  this  configuration  shown  in  Figure  5.3.  Separation  of  the  target  from 
clutter  is  entirely  based  on  the  difference  in  their  doppler  shifts.  Eventhough 
transient  response  is  no  longer  a  problem,  this  processor  suffers  from  several 
drawbacks.  Since  there  is  no  clutter  attenuation,  the  word  length  in  the  FFT 
processor  must  be  large  enough  to  accomodate  the  large  ground  clutter.  The 
absence  of  the  low  pass  filter  dictates  that  the  sampling  rate  at  the  FFT  processor 
be  maintained  at  5  Mil z .  The  doppler  bins  are  still  separated  by  500  Hz  and  hence 
only  100  bins  are  required  to  cover  the  doppler  range  of  interest.  This  means  that 
the  input  to  the  FFT  contains  10000  samples  (2  m  sec  at  5  MHz  sampling  rate)  while 
only  the  first  100  output  samples  are  required.  Thus,  alternate  and  more  efficient 
implementation  of  the  FFT  must  be  considered.  The  weighting  requirement  before  FFT 
is  even  more  critical  than  before  because  of  the  large  amount  of  clutter  entering 
the  processor.  The  high  sampling  rate  and  the  large  word  length  requirement  makes 
this  processor  very  unattractive. 

5.2.4  Configuration  D-4 

In  an  effort  to  reduce  the  word  length  requirement  at  the  FFT  processor,  a  delay 
line  canceller  is  introduced  in  this  configuration  shown  in  Figure  5.4.  Since  the 
clutter  spectrum  has  components  at  multiples  of  code  repetition  frequency,  the 
delay  in  the  canceller  must  be  an  Integer  multiple  of  the  code  period.  It  is 
shown  in  Figure  5.5  that  a  delay  of  two  code  periods  is  desirable  since  it  has 
better  gain  characteristics  for  the  doppler  frequencies  of  interest.  It  is  shown 
in  Appendix  B  that  the  transient  response  of  a  two  pulse  canceller  can  be  fully 
suppressed  by  discarding  the  first  and  the  last  output  pulses.  This  corresponds  to 
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Figure  5.4  Configuration  D-4 


Magnitude  Responses  of  Two  Delay  Line  Cancellers 


two  code  periods  in  the  beginning  and  at  the  end  of  the  output,  which  is  a  small 
portion  of  the  total  look  time.  Furthermore,  it  is  shown  in  Appendix  C  that 
eventhough  the  1TTI  filter  gain  is  not  the  same  for  all  dopplers,  the  signal-to- 
noise  ratios  in  the  various  FFT  doppler  bins  are  nearly  the  same.  However,  the 
signal  level  varies  from  one  bin  to  the  other  and  a  different  threshold  must  be 
chosen  for  each  bin.  The  major  disadvantage  of  this  configuration  is  the  fact 
that  the  FFT  processor  must  operate  with  an  input  sampling  rate  of  5  MHz. 

5.2.5  Configuration  D-5 

Due  to  the  presence  of  enormous  amounts  of  clutter,  large  word  lengths  are  required 
in  the  processor  until  the  clutter  rejecticn  filter.  It  is  therefore  advantageous 
to  perform  clutter  suppression  as  early  in  the  processor  as  possible.  In  this 
configuration  shown  in  Figure  5.6,  the  delay  line  canceller  is  implemented 
immediately  after  A/D  conversion.  It  is  shown  in  Appendix  D  that  the  two  pulse 
canceller  before  or  after  decoding  are  equivalent  provided  the  delay  is  an  integer 
multiple  of  the  code  period.  In  this  case,  the  performance  of  D-4  and  D-5  are 
identical.  However,  only  one  canceller  is  needed  in  D-5  as  opposed  to  one  in  each 
range  channel  in  D-4.  Also,  the  early  reduction  in  dynamic  range  could  lead  to  a 
more  economical  implementation. 

5.2.6  Configuration  D-6 

Shown  in  Figure  5.7,  this  configuration  is  the  same  as  D-3  except  that  a  low  pass 
filter  is  introduced  prior  to  FFT  processing  in  order  to  reduce  the  sampling  rate. 

Eventhough  the  number  of  samples  into  the  FFT  processor  is  greatly  reduced,  there 

I 

1 

is  no  loss  in  the  processor  gain  since  low  pass  filtering  provides  bandwidth  ] 

reduction  and  hence  an  increase  in  the  signal-to-noise  ratio.  Due  to  the  absence  j 

I 

of  a  clutter  rejection  filter,  large  amounts  of  clutter  enter  the  FFT  processor 
thus  requiring  a  large  word  length.  Furthermore,  due  to  the  sidelobcs  inherent 
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Configuration  D-5 


Figure  5.7  Configuration  D-6 


in  FFT  processing,  the  large  clutter  component  will  leak  into  other  doppler  bins 
unless  the  input  is  properly  weighted. 

5.2.7  Configuration  D-7 

This  processor,  shown  in  Figure  5-8,  combines  the  advantages  of  configurations 
D-4  and  D-6  in  that  both  a  clutter  rejection  filter  and  a  sampling  rate  reduction 
filter  are  utilized.  Here  again,  the  low  pass  filter  may  be  recursive  or  non¬ 
recursive,  and,  if  possible,  filtering  and  sampling  rate  reduction  may  be  combined 
by  implementing  a  fixed  window  transversal  filter.  Such  a  filter  contains  only 
one  multiplier  and  one  accumulator  in  each  range  channel  as  opposed  to  several 
multipliers  for  any  other  type  of  filter  realization.  As  in  the  case  of 
Configuration  D-2,  weighting  the  input  signal  is  required  to  reduce  interference 
between  two  target  returns  in  the  same  range  cell. 

5.2.8  Configuration  D-8 

This  is  exactly  the  same  as  the  previous  processor  except  that,  as  shown  in 
Figure  5.9,  the  clutter  rejection  filter  is  implemented  before  decoding.  Only 
one  delay  line  canceller  is  needed  as  opposed  to  63  in  Configuration  D-7. 

5.2.9  Conf i_g u ration  D-9 

The  main  motive  behind  this  configuration  is  the  possibility  of  using  a  very  small 
word  length  (as  small  as  one  bit)  at  the  decoder.  In  order  to  be  able  to  achieve 
this,  the  doppler  information  must  be  extracted  before  the  decoding  operation.  In 
this  configuration  shown  in  Figure  5.10,  the  incoming  signal  is  first  divided  into 
several  doppler  channels  and  range  resolution  is  later  achieved  in  each  doppler 
channel.  The  block  named  compensator  is  a  complex  weighting  of  the  signal,  and  it 
is  different  for  each  doppler  channel  as  described  in  Appendix  K.  Also  shown  in 
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Appendix  E  is  the  fact  that  this  configuration  is  identical  to  D-5  in  its 
performance.  The  gain  in  signal- to-noise  ratio  is  partially  achieved  by  the 
FFT  processor  and  the  rest  of  the  gain  is  attained  during  integration  following 
the  decoder. 

5.3  PROCESSORS  USING  PULSE  COMPRESSION 

The  pulse  compression  approach  is  an  attempt  to  convert  the  CW  signal  into  a 
pulsed  signal,  so  that  standard  pulse  radar  processing  techniques  may  be  used. 

The  pulse  compression  filter  is  a  matched  filter,  which  is  matched  to  one  period 
of  the  code.  In  the  absence  of  any  doppler,  the  output  of  the  matched  filter 
due  to  the  target  return  is  the  same  as  the  autocorrelation  of  the  code.  In  the 
presence  of  doppler  shift,  the  doppler  envelope  is  impressed  on  the  output  as 
shown  in  Figure  5.11.  However,  because  of  the  doppler  mismatch  (matched  filter 
is  matched  to  zero  doppler)  the  output  is  no  longer  the  autocorrelation  but 
features  larger  range  sidelobes,  thereby  reducing  range  resolution.  The.  mismatch 
also  reduces  the  main  peak  resulting  in  a  weaker  detection  performance.  Two 
configurations  PC-1  and  PC-2  are  shown  in  Figures  5.12  and  5.13,  the  only 
difference  between  the  two  being  the  location  of  the  clutter  rejection  filter. 

It  is  shown  in  Appendix  F  that  the  output  of  the  pulse  compression  processor  is 
not  the  same  as  that  using  the  decoder  except  at  zero  doppler.  It  must  be  noted 
that  the  implementation  of  the  processors  using  pulse  compression  is  considerably 
simpler  than  those  using  decoders  because  fewer  components  are  required.  It  must 
be  noted  too  that  a  waveform  ideally  suited  for  the  decoder  approach  may  not  be 
suitable  for  the  pulse  compression  approach  and  vice  versa. 
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12  Configuration  PC-! 


Configuration  PC-! 


5.4  HYBRID  PROCESSORS 


The  principal  limitations  of  the  digital  processors  are  the  conflicting  require¬ 
ments  of  high  sampling  rate  and  large  word  length,  while  their  advantages  are 
flexibility  and  easy  realization  of  coherent  integration  using  FFT  processors. 

The  hybrid  processors  are  designed  to  maintain  these  advantages  while  overcoming 
some  of  the  drawbacks  of  digital  processing  by  performing  a  portion  of  the 
processing  before  A/D  conversion.  Coherent  integration  capability  of  digital 
processing  is  utilized  in  Configuration  H-l  shown  in  Figure  5.14.  'Because  of 
clutter  attenuation  and  bandwidth  limiting  by  the  analog  filter,  both  the 
requirements  of  sampling  rate  and  word  length  are  made  less  stringent.  The 
analog  filtering  may  be  performed  either  at  video  on  the  in-phase  and  quadrature 
signals  or  at  IF.  If  IF  filtering  is  used,  the  in-phase  and  quadrature  components 
must  be  extracted  prior  to  A/D  conversion.  The  main  disadvantage  of  this  processor 
is  the  same  as  that  of  the  analog  processor  described  earlier,  that  is,  the  loss  in 
achievable  processor  gain  due  to  the  large  settling  time  of  the  notch  filter. 

This  problem  can  be  alleviated  by  using  Configuration  H-2,  shown  in  Figure  5.15, 
where  clutter  rejection  is  performed  in  the  digital  portion.  This  is  accomplished 
by  the  use  of  a  delay  line  canceller  in  which  the  transient  response  problems  arc 
completely  avoided.  However,  since  no  clutter  attenuation  takes  place  before  A/D 
conversion,  large  word  lengths  are  required  to  accomodate  the  expected  dynamic 
range.  This  problem  is  not  as  severe  as  in  the  all-digital  processors  because  of  the 
availability  of  A/D  converters  with  large  word  lengths  at  lower  sampling  rates. 

It  must  be  noted  here  that  eventhough  the  sampling  rate  has  been  reduced,  the 
number  of  A/D  converters  has  been  increased  from  one  to  one  in  each  range  channel. 
Superiority  of  hybrid  processing  is  accompanied  by  a  loss  in  flexibility  due  to 
the  fact  that  the  analog  portion  of  the  processor  cannot  be  easily  adapted  to  any 
changes  in  the  waveform  or  other  parameters. 


Figure  5.14  Configuration  H-l 


Figure  5.15  Configuration  H-2 


5.5  CONCLUDING  REMARKS 


The  principal  attractions  of  digital  processors  are  the  flexibility  they  offer  and 
the  possibility  of  obtaining  filter  characterics  which  could  not  be  realized  by 
analog  systems.  In  this  chapter,  several  configurations  were  proposed  as  possible 
candidates  for  digital  processors  in  a  phase  coded  CW  radar.  Each  processor  is 
based  on  one  of  two  different  concepts:  decoding  and  pulse  compression.  Even- 
though  several  of  the  configurations  are  mathematically  equivalent,  their  detection 
performance  may  vary  because  of  the  nonlinearity  inherent  in  quantization. 

Therefore,  each  of  th  processors  must  be  analyzed  in  detail  to  minimize  the 
cost  for  specified  performance  levels. 

Figure  5.16  shows  the  uncertainty  function  for  a  PN  sequence  of  length  63,  a  clock 
rate  of  5  MHz  and  a  look  time  of  2  m  sec.  As  can  be  seen  the  signal  exhibits  the 
desirable  properties  of  good  range  and  doppler  resolution  and  small  range  side  lobes. 
These  properties  can  only  be  attained  if  the  processor  consists  of  a  bank  of  filters 
each  matched  to  a  different  doppler  shift.  Configuration  D-3  is,  in  fact,  such  a 
processor  while  the  other  decoder  configurations  using  coherent  integration  are 
close  approximations.  However,  the  processors  using  pulse  compression  produce 
large  range  side  lobes  because  of  the  doppler  mismatch  in  the  pulse  compression 
filter  which  is  matched  to  zero  doppler.  This  can  be  seen  from  Figure  5.17  which 
shows  the  magnitudes  at  the  output  of  the  pulse  compression  filter  for  various 
dopplers.  The  input  to  the  pulse  compression  filter  is  a  doppler  shifted  phase 
code  at  zero  range  delay.  At  zero  doppler,  the  range  side  lobes  are  exactly  the 
same  as  in  the  decoder  approach,  since  the  pulse  compression  filter  is  exactly 
matched  to  the  input  waveform.  However,  as  the  doppler  shift  of  the  input  increases, 
the  mismatch  produced  manifests  itself  in  large  range  side  lobes  and  a  reduced  main 
lobe.  While  this  effect  is  not  pronounced  at  small  doppler  shifts,  the  degradation 
is  unacceptably  large  at  higher  doppler  frequencies.  It  must  also  be  noted  that  the 


5-23 


range  side  lobe  structure  depends  upon  the  initial  shift  in  the  code  (initial 
contents  of  the  shift  register)  unlike  in  the  decoder  approach.  The  severity 
of  the  range  side  lobes  is  entirely  dependent  upon  the  modulating  waveform 
and  it  is  necessary  to  analyze  several  waveforms  in  order  to  select  one  which 
is  suited  for  the  pulse  compression  processor. 

The  main  factor  that  restricts  the  performance  of  a  digital  processor  is  the 
effect  of  finite  word  lengths.  This  effect  can  be  reduced  by  expanding  the 
word  length,  but  an  upper  limit  on  the  word  length  is  imposed  by  the  high 
sampling  rate  necessary  to  avoid  errors  due  to  aliasing.  In  an  effort  to 
reduce  the  sampling  rate  required,  two  hybrid  processing  schemes  were  proposed 
where  some  flexbility  is  sacrificed  by  performing  a  portion  of  the  processing 
before  analog- to-digital  conversion. 


SECTION  6  -  SIMULATION  AND  RESULTS 


6.1  INTRODUCTION 

In  order  to  study  the  performances  of  the  various  configurations,  digital 
computer  simulations  of  the  processors  have  been  developed.  To  facilitate  a 

Monte  Carlo  analysis  of  the  performance,  computer  programs  have  also  been 

developed  to  generate  samples  of  synthetic  video  return  including  ground  clutter, 
noise  and  return  from  a  single  moving  target.  All  the  programs  are  written  in 
FORTRAN  compatible  with  the  CDC  6600  computer  system.  Program  listings  are 
included  in  Appendix  H. 

6.2  SYNTHETIC  VIDEO  RETURN 

Samples  of  the  video  return  are  normalized  so  that  the  power  in  the  target  return 

is  unity.  It  is  assumed  that  the  sampling  rate  is  equal  to  the  clock  rate  of  the 

PN  sequence.  The  total  return  includes  contributions  due  to  the  target,  noise 
and  clutter. 

6.2.1  Target  Return 

The  target  is  assumed  to  be  a  moving  target  with  a  dopplcr  frequency  f^  and  a 
range  delay  corresponding  to  k  clock  periods.  ihus,  if  the  target  range  is  R  and 
the  range  resolution  is  AR,  th.cn  the  range  delay  in  clock  periods  is 

k  =  Int  (R/AR)  modulo  (N)  (6.1) 

where  Int  (•)  denotes  "integer  part  of"  and  N  is  the  length  of  the  code.  If  the 
set  {p  ( i )  }  denotes  samples  of  tiie  periodic  code,  then  the  i1-*'  sample  of  the  in- 
phase  and  quadrature  components  of  the  target  return  arc  given  by 


xj(i)  =  p(i-k)  cos  (27rfdi6  +  '^t)  (6.2) 

Xq(i)  =  p(i-k)  sin  (2iTfdi6  +  i|>t)  (6.3) 

where  ^  is  a  random  phase  which  is  constant  over  a  look  period.  As  can  easily 
be  seen,  the  total  power  in  the  target  return  is  unity. 

6.2.2  Noise 

If  the  required  signal-to-noise  ratio  in  dB  is  SNR,  then  the  total  noise  power 
is 

a2  «  10-0.1  SNR  (6.4) 

This  noise  power  is  divided  equally  between  the  in-phase  and  quadrature  channels. 
The  noise  samples  in  each  channel  are  assumed  to  have  a  Gaussian  distribution  with 
zero  mean  and  a  variance  of  o2/2.  The  samples  are  independent  of  one  another  and 
the  noise  sequences  in  the  two  charnels  are  assumed  to  be  independent  of  each 

other. 


6.2.3  Ground  Clutter  Return 

If  the  required  signal-to-clutter  ratio  in  dB  is  SCR,  then  the  clutter  power  is 
A2  =  10"0-1  SCR  .  (6.5) 


Ground  clutter  return  is  assumed  to  have  a  steady  or  DC  component  and  a  random 
fluctuating  or  AC  component.  The  sample  of  clutter  return  in  the  in-phase 
and  quadrature  channels  are  generated  according  to  the  following  equations. 


cj(i)  =  A  p(i-k)  [  gQ  cos  i|>t.  +  ga  f1(i)]  (6.6) 

cg(i)  =  A  p(i-k)  [  gG  sin  i pc  +  ga  fq(i)l  (6.7) 
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where  { p ( ±) }  is  the  periodic  code  sequence, 


k  is  the  range  delay  of  the  clutter  source  in  units  of  clock  period. 


80 


_  ^1  ;  gQ  is  the  DC  component  of  the  clutter  power, 

'  W 


m ^  is  the  clutter  DC-to-AC  power  ratio. 


J — L_  ; 

*  2 ( 1+ni ) 


2 

2ga  is  the  AC  component  of  the  clutter  power. 


is  a  random  phase  angle  which  is  constant  over  a  look  period, 

and  (fj(i)}  and  {fg(i)}  are  samples  of  the  AC  clutter  component.  The  random 
sequences  {f-(i)l  and  {fg(i)l  are  assumed  to  have  a  Gaussian  distribution  with 
zero  mean  and  unit  variance  and  the  two  sequences  are  independent  of  each  other. 
However,  each  sequence  is  assumed  to  be  highly  correlated  and  to  have  a  Gaussian 
spectrum  with  zero  mean  and  a  standard  deviation  of  Hz.  Samples  of  the 
random  clutter  component  are  generated  by  introducing  the  desired  correlation  to 
an  independent  random  sequence.  This  is  accomplished  using  a  linear  system,  and 
the  process  is  explained  in  Appendix  G. 


Various  types  of  clutter  may  be  produced  by  proper  use  of  the  above  procedure. 

Pure  DC  clutter  resulting  from  static  reflectors  such  as  buildings  crn  be  obtained 
by  setting  ga  to  zero  and  g0  to  one.  Antenna  spillover  can  be  simulated  as  a  DC 
clutter  with  zero  range  delay,  that  is,  k=0.  Distributed  clutter  is  generated  by 
adding  the  clutter  returns  from  several  range  cells.  The  return  from  each  range 
cell  must  take  into  account  the  power  variation  according  to  the  inverse  of  the 
fourth  power  of  range,  and  the  delay  corresponding  to  the  particular  cell. 


6.3  PROCESSOR  SIMULATION 


Since  the  proposed  digital  processors  contain  similar  components  but  in  different 
orders,  each  of  the  blocks  was  separately  simulated.  This  includes  A/D  converter, 
decoder,  pulse  compression  filter,  delay  line  canceller,  weighting  function,  digital 
filter  for  clutter  rejection  and  sampling  rate  reduction,  FFT  processor,  and  non¬ 
coherent  integrator.  Simulation  of  the  digital  processors  is  not  complete  in  that 
only  the  quantization  errors  due  to  A/D  conversion  can  be  simulated.  Finite  word 
length  effects  elsewhere  in  the  processor  including  the  fast  Fourier  transformation 
have  not  been  included.  The  notch  and/or  sampling  rate  reduction  filter  is  assumed 
to  be  a  recursive  filter  and  is  implemented  as  a  cascade  of  second  and  first  order 
sections  to  minimize  the  effects  of  coefficient  quantization.  Each  of  the  first  and 
second  order  sections  is  realized  in  the  canonical  structure.  It  should  be  noted 
that  processor  Configuration  D-9  cannot  be  simulated  with  the  programs  that  have  been 
developed  because  of  its  many  differences  in  comparision  with  the  other  configurations. 

6.4  A  SIMULATION  EXPERIMENT 

The  simulations  developed  were  exercised  for  a  typical  set  of  parameters.  The  object 
of  the  experiment  was  twofold:  (i)  to  study  the  feasibility  of  digital  processing 
when  the  target  return  is  much  smaller  than  the  quantization  interval  of  the  A/D 
converter  and  (ii)  to  compare  the  performance  of  a  digital  processor  with  that  of 
the  analog  processor  described  in  Chapter  3. 

Configuration  D-8  was  chosen  as  the  candidate  digital  processor.  Figure  6.1  shows  the 
magnitude  response  of  a  fourth  order  Chebyshev  low  pass  filter  used  for  sampling  rate 
reduction.  This  filter  lias  5  dB  ripples  in  its  pans  band  but,  using  the  same 
argument  as  with  the  delay  line  canceller,  this  should  not  affect  the  output  signnl- 
to-noise  ratio  due  to  the  narrow  doppler  filters  provided  by  FFT  processing.  The 
output  of  this  low  pass  filter  was  resampled  at  80  KMz  as  opposed  to  64  KHz  suggested 


in  Figure  5.9.  The  input  samples  to  the  FFT  processor  were  weighted  with  a 
Hamming  function  shown  in  Figure  6.2.  The  simulation  of  the  analog  processor 
is  exactly  similar  to  Configuration  D-l  with  one  difference.  The  non-coherent 
integration  is  performed  at  the  original  sampling  rate  instead  of  a  reduced 
sampling  rate  to  provide  a  better  simulation  of  the  analog  integration. 

The  analog  notch/low  pass  filter  was  modelled  by  a  10th  order  Butterworth  low 
pass  filter  in  cascade  with  an  8th  order  Butterworth  high  pass  filter.  The 
magnitude  response,  shown  in  Figure  6.3,  was  found  to  be  a  satisfactory  approximation 
to  the  response  of  the  analog  filter.  Assuming  a  look  time  of  2  milliseconds,  a 
non-symmetrical  Taylor  weighting  function  (with  n  =  10  and  -50  dB  sidelobes),  shown 
in  Figure  6.4,  was  applied  to  the  input  prior  to  filtering.  As  explained  in 
Chapter  3,  the  transient  response  must  be  allowed  to  settle  down  before  starting 
integration.  In  this  experiment,  the  non-coherent  integration  was  started  after 
1.33  milliseconds,  i.e.,  out  of  the  2  milliseconds  look  time,  integration  was 
performed  only  over  0.67  millisecond. 

Table  6.1  shows  the  radar  parameters  used  in  the  simulation.  For  these  set  of 
typical  radar  parameters,  the  signal- to-noise  ratio  at  the  input  to  the  receiver 
can  be  computed  to  be  -20  dB  using  the  radar  range  equation.  A  complete  list  of 
other  simulation  parameters  is  given  in  Table  6.2. 

6.5  SIMULATION  RESULTS 

The  outputs  in  range  channels  32,  28  and  36  were  computed  for  both  the  analog  and 
the  digital  processor  using  the  parameters  given  above.  These  channels  represent 
the  in-range  and  two  out-of-range  channels  all  of  which  contain  clutter  and  noise. 

A  Monte-Carlo  analysis  of  30  samples  was  performed  to  determine  the  means  and 
variances  of  the  outputs.  These  results  arc  given  in  Table  6.3  where  the  output 


6-6 


•But  te  rwor  th  iFi  I  tec Magrjit  ude  jRe  sports  e  | 


Table  6.1  Radar  Parameters  Used  in 

Transmitted  Power 
Transmit  Antenna  Gain 
Receive  Antenna  Gain 
Wavelength  (X-band) 

Target  Cross  Section  (Steady) 

Target  Range 
Length  of  PN  Sequence 
Clock  Rate 
Code  Period 
Target  Range  Rate 
Doppler  Frequency 
Receiver  Bandwidth 
Noise  Figure 


Simulation 

10  watts 
38  dB 
40  dB 

0.03  meter 
1  sq .  meter 
15000  meters 
63 

5  MHz 

12.6  psec. 

240  meters/sec. 
16  KHz 
5  MHz 
5  dB 


Losses 


9.5  dB 


Table  6.2  Simulation  Parameters 

Target:  1  sq.  meter  steady  target  in  range  cell  32 

In-range  Clutter:  Signal-to-clutter  ratio  =  -70  dB 

DC-to-AC  power  ratio  =  0.8 

Out-of-range  Clutter:  Signal-to-clutter  ratio  =  -70  dB 

DC-to-AC  power  ratio  =0.8 
Appears  in  range  cell  28 

Fixed  Clutter:  Signal-to-clutter  ratio  =  -60  dB 

Appears  in  range  cell  32 

Clutter  Parameters:  Wooded  Terrain 

Wind  Velocity  =  18  knots 
Clutter  doppler  spread  a-p  =  8  Hz 

Noise:  Signal- to-noise  ratio  =  -20  dB 

Sampling  Rate:  5  MHz 

A/D  Converter  Input  Range:  -18000  to  18000  volts 

Look  Time:  2  milliseconds  for  analog  processor 

1.65  milliseconds  for  digital  processor 

Integration:  Square  law  detection  and  integration  for  0.67 

millisecond  in  analog  processor 

128  sample  FFT  in  digital  processor  followed 
by  a  square  law  detector 
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Table  6.3  Simulation  Results 


Range  Bin 

28 

Range  Bin 

32 

Range  Bin 

36 

Analog  Processor 

Mean 

8692.04 

17612.51 

8472.89 

Variance 

4 . 6585xl06 

6.5340xl06 

1.4319xl06 

Digital  Processor 

36  Bits 

Mean 

0.5044 

33.0759 

0.4370 

Variance 

0.1752 

29.5580 

0.1726 

Digital  Processor 

10  Bits 

Mean 

1.2152 

32.4057 

0.9593 

Variance 

2.9316 

91.9885 

0.5703 

Digital  Processor 

13  Bits 

Mean 

0.5493 

33.3664 

0.4443 

Variance 

0.1910 

29.4330 

0.2020 
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of  the  digital  processor  corresponds  to  the  output  in  the  doppler  channel  in 
which  the  target  is  present.  The  digital  processor  was  analyzed  for  three 
different  A/D  converter  word  lengths:  36  bits,  10  bits  and  13  bits.  The  36 
bit  word  length  was  included  for  validation  purposes  since  the  quantization 
effects  at  this  word  length  are  inconsequential. 

6.6  ANALYSIS  OF  THE  RESULTS 

In  the  case  of  the  digital  processors,  the  mean  of  the  output  is  the  variance  before 
detection  since  a  square  law  detector  is  being  used.  Since  the  original  bandwidth 
is  5  MHz  and  the  detection  bandwidth  is  620  Hz  (based  on  the  look  time),  the  expected 
gain  in  signal-to-noise  ratio  is  39  dB.  Thus  the  expected  signal-to-noise  ratio  at 
the  output  without  quantization  effects  is  19  dB.  Using  the  results  for  36  bits  and 
range  bins  28  and  32,  it  can  be  show  that  the  signal- to-interference  ratio  is  18  dB. 
This  agrees  with  the  predicted  value  and  suggests  that  the  large  amount  of  input 
clutter  has  an  insignificant  effect  on  the  output. 

With  a  13  bit  A/D  converter,  the  peak  signal  amplitude  is  four  times  smaller  than 
the  quantization  interval  corresponding  to  a  signal-to-quantization  error  ratio  of 
-5  dB.  The  effect  of  this  is  to  decrease  the  input  signal-to-noise  ratio  from  -20  dB 
to  -20.13  dB .  The  results  show  that  the  output  signal-to-noise  ratio  is  17.75  dB, 
which  is  a  fraction  of  a  dB  below  that  for  the  case  where  the  word  length  is  36  bits. 
It  should  be  noted  that  with  iB  bits,  the  standard  deviation  of  the  input  noise  is 
about  twice  the  quantization  interval. 

When  the  word  length  is  reduced  to  10  bits,  the  input  noise  standard  deviation  is 
one-fifth  the  quantization  interval.  The  corresponding  signal-to-quantization 
error  is  -23.15  dB  which  has  the  effect  of  increasing  the  input  signal-to-noise 
ratio  from  -20  dl!  to  -24.7  dB.  From  the  results  it  can  be  shown  that  the  output 
signal-to-noise  ratio  is  14  dB  which  is  4  dB  less  than  the  unquantized  case. 
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Due  to  the  noncoherent  integration  and  the  fact  that  clutter  residues  persist 
despite  the  weighting  and  delaying  the  start  of  integration,  it  is  not 
straightforward  to  interpret  the  signal-to-interference  ratio  at  the  output  of 
the  analog  processor.  The  most  desirable  way  to  analyze  this  processor  is  to 
perform  an  extensive  Monte  Carlo  analysis  to  ascertain  the  probabilities  of 
detection  and  false  alarm.  Such  a  simulation  requires  an  enormous  amount  of 
computation.  For  purposes  of  this  preliminary  comparison,  probabilities  of 
detection  and  false  alarm  were  computed  by  using  approximate  probability  density 
functions.  Probability  densities  with  and  without  a  target  were  approximated  by 
a  non-central  and  a  central  chi-squared  distribution  respectively.  The 
appropriate  parameters  of  these  distributions  were  derived  using  the  results  of 
the  simulation.  Then,  by  varying  the  threshold,  a  curve  of  probability  of  false 
alarm  versus  probability  of  detection  was  developed  for  each  processor.  This 
set  of  curves  is  shown  in  Figure  6.5. 

6.6  CONCLUSIONS 

The  results  of  the  simulation  show  that  it  is  indeed  possible  to  extract  the 
signal  using  a  digital  processor  whose  word  length  is  such  that  the  signal  is 
considerably  smaller  than  a  quantization  interval.  Furthermore,  detection  is 
possible  even  when  the  input  noise  standard  deviation  is  smaller  than  the 
quantization  interval.  Despite  the  large  transients  caused  by  clutter,  i t  is 
shown  that  the  analog  processor  can  be  used  if  a  proper  weighting  function  is 
used.  However,  for  the  set  of  parameters  used,  a  look  time  of  2  milliseconds 
docs  not  seem  to  provide  an  adequate  detection  performance  by  the  analog 
processor.  Finally,  for  a  given  look  time,  the  digital  processor  gives  a 
considerably  better  detection  performance.  This  is  attributable  to  the  absence 
of  undcsriablo  transients  and  the  use  of  coherent  integration. 
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SIX  T I  ON  7  -  SUMMARY,  CONCLUSIONS,  AND  FUTURE  EFFORTS 

The  principal  problem  that  must  be  addressed  in  designing  a  signal  processor  for 
a  phase  coded  CW  radar  is  the  large  dynamic  range  requirement  that  is  caused  by 
antenna  spillover  and  close-in  clutter.  Since  the  clutter  return  could  be  60-80  dB 
over  the  target  return,  the  clutter  rejection  filter  must  be  designed  to  have  a 
very  short  settling  time.  Otherwise,  the  transient  response  due  to  the  large 
clutter  return  will  have  the  effect  of  reducing  the  signal- to-interference  ratio 
thereby  deteriorating  the  detection  performance  of  the  radar.  Since  digital 
clutter  rejection  filters,  such  as  delay  line  cancellers,  can  be  designed  to  have 
very  short  rettling  times,  the  possibility  of  using  digital  processing  must  be 
given  serious  consideration.  Another  attractive  feature  of  digital  processors  is 
that  they  are  considerably  more  flexible  than  analog  processors.  The  main 
disadvantage  of  digital  processing  is  the  introduction  of  quantization  errors 
due  to  the  finite  word  length. 

For  a  given  dynamic  range  requirement,  quantization  errors  can  be  reduced  by 
increasing  the  word  length.  But  the  high  sampling  rate  made  necessary  by  the  use. 
of  a  wideband  signal  imposes  a  limit  on  the  word  length  that  can  be  obtained  from 
currently  available  analog-to-digital  converters.  At  the  dynamic  range  and 
sampling  rate  of  interest  in  this  application,  the  largest  word  length  available 
has  a  quantization  interval  which  is  much  bigger  than  the  amplitude  of  the  target 
return.  This  problem  can  be  partially  alleviated  by  using  hybrid  processors 
where  some  flexibility  of  digital  processing  is  sacrificed  to  relax  the  stringent 
and  conflicting  requirements  of  large  dynamic  range  and  high  sampling  rate. 

Several  analog,  digital  and  hybrid  configurations  were  proposed  in  this  report 
as  candidates  for  signal  processors  in  a  binary  phase  coded  CW  radar.  Synthetic. 
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input  signals  including  returns  due  to  clutter,  target  and  noise  were  used  along 
with  digital  computer  simulations  of  the  processors  to  show  the  feasibility  of 
digital  processing.  It  was  shown  that  is  was  possible  to  detect  a  target  when 
the  target  return  is  much  smaller  Lhan  the  quantisation  interval  of  the  A/D 
converter.  The  target  was  detectable  even  when  the  input  noise  was  smaller  than 
the  quantization  interval.  It  was  also  shown  that  the  analog  processor  could  be 
used  if  the  transient  response  of  the  clutter  rejection  filter  were  controlled  by 
modifying  the  input  signal  by  a  suitable  weighting  function. 

It  must  be  remarked  that  this  study  is  incomplete  and  that  the  most  significant 
conclusion  is  that  digital  processing  is  feasible  despite  the  rigid  sampling  rate 
and  dynamic  range  requirements.  The  following  is  a  partial  list  of  topics  that 
must  be  investigated  before  selecting  the  optimal  processors. 

1.  Effects  of  arithmetic  roundoff  and  quantization  of  filter  coefficients 
should  be  studied.  Since  no  adequate  theoretical  tools  are  available, 
this  must  be  performed  through  computer  simulation. 

2.  Each  of  the  configurations  should  be  analyzed  to  determine  dynamic 
range  requirements  at  different  stages  of  the  processor.  This  is 
required  in  order  to  be  able  to  use  the  available  word  length  in  the 
most  efficient  manner. 

3.  The  simulation  programs  must  be  modified  for  use  in  the  SEE  computer 
and/or  the  AP  12015  array  processor.  This  will  facilitate  the 
computation  of  detection  and  false  alarm  probabilities  using  an 
exhaustive  Monte-Carlo  analysis. 

4.  The  feasibility  of  using  a  fixed  window  transveral  filter  for  sampling 
rate  reduction  must  be  investigated.  This  is  important  since  such  a 


filter  can  be  implemented  with  much  less  hardware  than  a  recursive 


filter . 

5.  Adaptive  thresholding  by  the  use  of  constant  false  alarm  rate  (CFAR) 
processors  should  be  examined.  This  technique  is  useful  in  avoiding 
false  alarms  in  the  presence  of  large  unexp'-e t ed  interference  sources. 

6.  MT1  filters  other  than  delay  line  cancellers  can  be  used  for  clutter 
rejection.  These  can  be  of  either  recursive  or  non- recurs i ve  variety. 

In  the  latter  case  proper  initialization  is  necessary  to  avoid  transient 
effects. 

7.  The  study  must  be  expanded  to  include  the  effects  of  bandlimiting  on 
the  code  and  the  effects  of  scanning  modulation  where  a  continuous 
scan  is  used  instead  of  a  step  scan. 

8.  The  processor  must  be  modified  to  operate  in  clutter  environments 
other  than  ground  clutter.  This  entails  the  use  of  notch  filters  for 
rejection  of  weather  clutter  and  chaff. 

9.  The  typical  s >t  of  parameters  used  in  this  study  yields  an  unambiguous 
range  of  1890  meters.  One  way  to  remove  this  ambiguity  is  to  increase 
the  length  of  the  I’N  sequence. 

This,  however,  reduces  the  code  repetition  frequency  which  in  turn 
decreases  the  separation  between  blind  speeds,  produced  by  the  clutter 
rejection  filter.  To  avoid  blind  speeds,  some  method  analogous  to 
staggered  FRF  must  he  employed.  The  extension  of  staggered  FRF  to  phase 
coded  CM  radars  is  not  obvious  and  needs  further  investigation.  The 
enlarged  length  of  the  sequence  also  increases  the  number  of  range  channels 
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and  hence  the  complexity  of  the  processor.  Therefore,  alternate 
schemes  for  unambiguous  range  measurement  must  be  examined. 

10.  The  maximum  length  binary  sequence  does  not  seem  to  be  ideally  suited 
for  the  pulse  compression  approach.  It  is  therefore  necessary  to 
study  other  types  of  waveform  to  find  one  which  provides  a 
satisfactory  performance  with  the  pulse  compression  processor.  It 
must  be  remembered  that  the  pulse  compression  processors  are 
considerably  easier  to  implement  than  those  using  decoders. 

11.  Even  though  some  of  the  digital  processor  configurations  are 
mathematically  equivalent,  their  performances  may  not  be  identical 
because  of  the  nonlinearity  involved  in  quantization.  It  is 
therefore  necessary  to  evaluate  all  the  configurations  to  select  a 
processor  which  offers  the  best  compromise  between  performance  and 
cost  of  implementation. 
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APPENDIX  A 


RESPONSE  OF  THE  FFT  PROCESSOR 

Let  the  input  to  the  processor  be  a  sampled  complex  exponential, 

Xjj  =  exp  ( 2tt  j  f  n  A) ,  n=0,  1,  .  .  ,  N-l  (A.l) 

where 

f  =  frequency 

N  =  number  of  input  samples, 

and  A  =  sampling  interval.' 

The  frequency  resolution  of  the  FFT  processor  is  the  reciprocal  of  the  input 
duration,  i.e., 

F  =  1/NA  ,  (A. 2) 


and  the  kth  output  sample  corresponds  to  a  frequency  kF.  The  frequency  response 
of  the  kth  frequency  bin  will  be  defined  as  the  magnitude  of  the  kth  output 
sample  as  a  function  of  the  input  frequency  f.  Tins  output  sample  is  given  by 


(O 


1 

N 


N-l 

E 


n=0 


(A.  3) 


where 

w  =  exp  (-2nj/N) 

Using  (A.l)  and  (A. 2),  (A. 3)  reduces  to 
,  N-l 

Xk  =  N  °Xp  (2l'JnAfr) 

n=0 

A-l 


(A.  A) 


(A. 3) 


where 


fr  =  f  -  k  F. 


(A.  6) 


The  summation  in  (A. 5)  is  a  finite  geometric  series  and  can  be  written  as 


xkCf)  -  *tr.)  -1 

K  N  exp  (2ttj  /.fr)  -1 


(A.  7) 


The  above  equation  can  be  rearranged  to  give 


xk(f),=  A 

k  N 


sin  ttN  Afr 
sin  n  A  fr 


(A. 8) 


For  large  values  of  N  and  small  values  of  f  ,  the  frequency  response  can  be 
approximated  as, 


xk(f) 


-  sa  (ir  NA  fr) 


(A.  9) 


where  the  function  sa(x)  denotes  sin  x/x. 


An  examination  of  the  above  expression  reveals  the  following  properties: 

(i)  the  response  of  the  kth  bin  is  centered  at  frequency  kF, 

(ii)  the  frequency  response  contains  large  sidelobes,  and 

(iii)  there  is  considerable  overlap  between  the  responses  of  neighboring 
bins.  In  fact,  at  the  frequency  where  the  two  responses  cross, 
the  gain  is  about  4  dB  below  the  maximum  gain. 

If  the  input  is  an  independent,  zero  mean  noise  sequence  with  variance  1,  it  can 
be  easily  shown  that  the  output  in  each  bin  of  the  FFT  processor  is  a  zero  mean 
noise  whose  variance  is  1/N.  Since  the  input  signal -to-noise  ratio  is  unity,  the 
gain  in  signal- to-noise  ratio  is  the  same  as  the  output  signal- to-noise  ratio 
which  is. 
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RESPONSE  OF  DELAY  LINE  CANCELLER 

The  input  is  assumed  to  be  a  finite  duration  signal, 

r(t)  =  x(t)  Wl(t)  (B-l) 

where  x(t)  is  an  infinite  duration  signal  and  w^(t)  is  a  window  function  given 
by 

w^t)  *»  1  f  r  0  <t<  T 

=  0  otherwise  (B.2) 

The  output  of  the  two  pulse  delay  line  canceller  is 

y(t)  =  x ( t)  -  x(t-A)  (B.3) 

and  the  impulse  response  of  the  filter  is 

h(t)  -  6 (t)  -  6(t-A)  .  (B.4) 

The  output  can  therefore  be  written  as 

y(t)  =  x ( t)  WjCt)  -  x(t-A)  w-^  (t-A)  (B.5) 

and  exists  in  the  interval  0<t<  T  +  A. 

To  avoid  the  effects  of  transients,  a  portion  of  this  signal  is  discarded. 
Specifically,  two  segments  of  y(t)  in  the  intervals  0<t<A  and  T<t<T  +  A  arc  not 


B-.1 


used  in  further  processing.  This  is  identical  to  multiplying  y(t)  by  a  window 
function  W2(t)  such  that, 


W2(t)  =  1  for  A<t<T 

=  0  otherwise  •  (B.6) 

But  from  the  definitions  of  wi(t)  and  W2(t),  it  can  be  easily  seen  that 

Wf(t)  W2(t)  =  W2(t),  (B.7) 

and  w^(t-A)  W2(t)  =  W2(t)  (B.8) 

Therefore  the  final  output  is 

z(t)  =  [x ( t)  -  x(t-A) ]  W2(t)  (B.9) 

Thus,  the  effect  of  discarding  the  transients  is  equivalent  to  passing  the 
infinite  duration  signal  through  the  canceller  followed  by  windowing.  The 
importance  of  this  is  in  regard  to  ground  clutter.  If  the  transients  are  not 
discarded,  the  effect  of  finite  observation  time  is  to  smear  the  spectrum  of  the 
input  signal.  Therefore,  the  low  frequency  ground  clutter  spectrum  expands  into 
the  passband  of  the  filter.  Discarding  the  transients  imnroves  clutter  rejection 
performance  since  spectrum  smearing  takes  place  after  the  filtering  which  is 
designed  to  attenuate  the  low  frequency  clutter  components  and,  in  particular,  to 
suppress  the  DC  clutter  completely.  It  must  be  noted  that  this  analysis  is  valid 
only  for  step  scan  systems  since  the  input  amplitude  has  been  assumed  constant. 
The  effect  of  discarding  transients  is  not  so  great  in  continuously  scanned 
systems  since  the  antenna  pattern  automatically  produces  a  spectral  widening  of 
the  input  signal. 
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APPENDIX  C 

COHERENT  INTEGRATION  FOLLOWING  MTI 


The  input  sampled  signal  is  assumed  to  be  of  the  form, 

r(n)  =  x(n)  +  p(n),  n  =  -1,  0,  1 . N-l  (C.l) 

where  x(n)  =  exp(27rj  f  nA)  (C  .2) 

and  p(n)  is  a  zero  mean,  independent  noise  sequence  with  unit  variance.  Notice 
that  the  input  signal-to-noise  ratio  is  unity.  It  is  assumed  that  this  signal 
is  passed  through  a  two  pulse  canceller  followed  by  an  FFT  processor.  The  signal 
component  of  the  output  of  the  two  pulse  canceller  after  discarding  transients  is, 

y (n)  =  x(n)  -x(n-l) 

=  exp (2-rr j  f  nA)  [l-exp(-2irj  fA)  ]  (C  .  3) 

If  this  sequence  is  passed  through  an  FFT  processor,  the  kth  output  sample  is 
i  N"1  nk 

Y(k)  =  I  S  y(n)w  (C,4) 

n=0 


where  w  =  exp(-2Trj/N)  . 

From  (C.3)  and  (C.4),  one  can  write 

N-l 

Y  (k)  =  [l-exp(-2irj  fA)]  —  exp(2Ti.i  f  nA)wn 

n=0 

N-l 

=  (l-exp(-27rjfA)  ]  —  ^2  exp(2wj  n  Afr  ) 

n=0 


f 


(C.5) 


where  f  =  f  -  kF 

and  F  =  1/NA  . 


It  follows  that. 


Y(k) 


2  simr  fA 


sin  it  N  A  f  r 
Nsin  7t  A  fr 


The  noise  samples  at  the  output  of  the  canceller  are 


m(n)  =  p(n)  -p(n-l),  n=0,  1,  ...  ,  N-l. 


Obviously,  this  is  a  zero  mean  noise  sequence  but  is  not  an 
due  to  the  correlation  introduced  by  the  canceller.  In  fact 
shown  tha  t 


E [m(n)m(i) ]  =  2 

=  -1 

=  0 


if  n  =  i 


if 


n-  i 


1 


otherwise 


where  E[.]  denotes  statistical  expectation. 


The  kth  output  noise  sample  of  the  FFT  processor  is, 
N-l 


M(k)  =  I  £  m(n)wnk  . 
n=0 


This  noise  sample  has  zero  mean  and  a  variance  given  by, 


E  [  |m  (k )|  2]  =  i  E 


N-l  N-l 


y  ^  y  ^  m(n)m(i  )wkn  w 
n=0  i-0 


C-2 


(C  .6) 

(C .  7) 

(C .  8) 

(C  .9) 

independent  sequence 
,  it  can  be  easily 


(C.10) 

(C.ll) 

(C .12) 


Using  (C.10),  the  above  equation  can  be  reduced  to 


E[  |M(k)|  2]  =  i  e 

N 


N-l 


m2 (n)-hn(n)m(n-l)w-k+m(n)m(n-l) 

n=0 


— 2  [2N  -  N(w^  +  w  ] 

N 


=  4  sin2(irk  FA)/N. 


(C.13) 


The  total  gain  in  signal-to-noise  ratio,  which  is  the  same  as  the  output 
signal-to-noise  ratio  is, 


4  sin^  tt f  A  sin  ttNa  <~r 

4  sin2  (likFA)  /N  Nsin  n  A  fr 

In  the  doppler  bin  containing  the  target,  f  -  kF  and  for  large  values  of  N, 

(C.14)  can  be  approximated  as, 

SNR0  ^  N  [sa(Tr  N  A  fr)  ]2  (C.15) 


2 

(C.14) 


which  is  the  same  as  that  obtained  by  coherent  integration  without  the  delay 
line  canceller.  Thus,  introducing  a  two  pulse  canceller  ahead  of  the  FFT 
processor  has  little  effect  on  the  improvement  in  signal-to-noise  ratio. 
However,  the  signal  levels  in  different  frequency  bins  are  not  the  same  and 
each  bin  must  be  assigned  a  different  threshold  which  takes  into  account  the 
frequency  response  of  the  delay  line  canceller. 


APPENDIX  D 


MTI  BEFORE  AND  AFTER  DECODING 


The  received  signal  consists  of  K  periods  of  the  code  and  there  are  N  samples  in 
each  code  period.  These  samples  are  denoted  as  r(i),  i=0,  1,  .  .  .  ,  KN-1 
Let  this  set  be  partitioned  into  K  segments  of  N  samples  (one  code  period)  each 
such  that, 

r^Cn)  =  r(kN  +  n)  for  k  =  0,  1,  .  .  .  ,  K-l, 

and  n  =  0,  1 . N-l.  (D.l) 

The  output  of  the  decoder  in  range  channel  corresponding  to  a  delay  of  m  samples 
is  obtained  by  multiplying  r (i )  by  a  periodic  repetition  of  Cm(n)  where  Cm(.) 
represents  a  period  of  the  code  shifted  by  m  samples.  Thus,  if  the  output  y(i) 
is  also  partitioned  as  above,  then 

yk(n)  =  rR(n)  Cm(n)  .  (D.2) 

If  this  sequence  is  then  passed  through  a  two  pulse  canceller  whose  delay  is  an 
integer  multiple  of  the  code  period,  i.e.,  the  delay  is  pk  samples,  the  output  of 
the  canceller  in  partitioned  form  is 

zk(n>  =  yk<n)  -  >'k-p(n) 

=  [rk(n)  -  rk_p(n) ]  Cm(n)  .  (D.3) 

If  the  MTJ  is  performed  before  decoding,  the  output  of  the  two  pulse  canceller 
in  partitioned  form  is 

xk (n)  =  rk(n)  -  rk_p  (n)  .  (D.4) 


D-l 


If  the  sequence  (n)  is  then  decoded,  the  decoder  output  in  the  range  channel 
corresponding  to  a  delay  of  m  samples  is, 

vk(n)  =  xk  (n)  '  cm(n) 

=  [rk(n)  -  rk-p  (n) ]  Cm( n)  .  (D.5) 

From  (!) .  0  and  (D.5)  it  is  seen  that  interchanging  the  MTI  and  decoding  operations 
has  no  effect  on  the  output.  However,  it  must  be  noted  that  this  result  is  based 
on  the  assumption  that  the  canceller  delay  is  an  integer  multiple  of  the  code 
period.  The  result  is  not  valid  for  other  delays  and  for  recursive  MTI  filters. 
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APPENDIX  E 


OUTPUT  OF  PROCESSORS  USING  DECODERS 

In  order  to  demonstrate  the  equivalence  of  Configurations  D-5  and  D-9, 
expressions  for  their  outputs  are  developed  in  this  Appendix.  In  the  development 
that  follows,  the  delay  line  canceller  is  ignored  since  its  effects  on  the  two 
configurations  are  identical.  The  input  signal  is  denoted  by  the  set  of  samples 
(r(i),  i  =  0,  1,  .  .  .  ,  KN-1  }  where  K  is  the  number  of  code  periods  in  the  look 
time  and  ’!  is  the  number  of  samples  in  each  code  period. 

In  Configuration  D-5,  the  output  of  the  range  channel  corresponding  to  a  delay  of 
m  samples  is  obtained  by  multiplying  the  input  by  a  periodic  repetition  of  Cm(n) 
where  Cm(.)  represents  a  period  of  the  code  shifted  by  m  samples .  Thus,  the 
decoder  output  is  given  by, 

K-l 

ym(i)  =  r(i)  ^2  CmCi-kN)  (E.l) 

k=0 

The  nth  outpuL  sample  following  an  NK  sample  Fl'T  is 
NK-1 

Ym(n)  =  S  ym(1)"'Jn  (E-2) 

i=0 

where  w  =  exp(-  2 ’■  j /’NK)  .  (E.l) 

Combining  (E.l)  and  (E.2), 

NK-1  K-l 

Yn>>  =  J2  lL  r(i)  cm(1-kN)win 
i=0  k=0 


E-  1 


I 
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Each  of  the  K  sequences  (U(p,n),  p  =  0,  1,  .  .  .  , 


N-l}  is  doppler  compensated, 


X(p,n)  =  U(p,n)  wnp  (E.10) 

The  output  of  the  processor  in  range  cell  in  and  doppler  cell  n  is  then  computed 
by  decoding  followed  by  summation, 

N-l 

Vm(n)  =  £  X(p,n)  Cm(p)  (E.  11) 

P=0 

N-l 

=  £  Zp(n)  wnp  (^(p) 

p=0 

N-l  K-l 

=  J  J  zP<k>wl  wllP  cm(P> 

p~0  k=0 

N-l  K-l  , M 

-  £  E  »"kN  »"p  q.(p> 

P=0  k =0 
N-l  K-l 

=  iC  2  r(kN+p)  cm(p)  wIlp  w"kN  '  (E.12) 

p=0  k=0 

Comparison  of  (E.5)  and  (E.12)  proves  the  hypothesis  that  the  outputs  of 
Configurations  (D-5)  and  (D-9)  are  identical. 
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APPENDIX  F' 


OUTPUT  OF  PROCESSORS  USING  PULSE  COMPRESSION 


Referring  to  Configuration  PC-2  and  ignoring  the  delay  line  canceller  as  in 
Appendix  Et  the  output  of  the  pulse  compression  filter  in  range  channel  m  is 
given  by 

N-l 

ym(k)  =  £  c^pMkN+p)  .  (F.i) 

p=0 

If  the  K  sample  FFT  oT  this  sequence  is  computed,  the  output  sample  in  the  nth 
bin  is 

K-i 

V")  -  E  y»(k)  “i  (F-2) 

k=0 

Here  w^  =  exp(-2Trj/K) 

=  w”  (F.3) 

where 

w  =  exp(-2irj/NK)  (F.4) 

Therefore, 

K-l  N-l 

Vn)  -  2  2  r(kN+P}  CB<P>W“kN  •  (F-5) 

k=0  p=0 

Comparison  of  (F.5)  and  (E.12)  reveals  the  similarity  between  the  decoder  and  the 
pulse  compression  approaches,  the  only  difference  being  the  term  due  to  doppler 
compensation. 
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APPENDIX  G 


CLUTTER  GENERATION 


It  is  desired  to  generate  samples  of  the  random  clutter  components  (fj(i)}  and 
{ f q(i) }  at  sampling  intervals. of  6.  It  is  assumed  that  the  sequence  can  be 
described  by  a  Gaussian  distribution  with  zero  mean  and  unit  variance.  These 
sequences  are  required  to  be  highly  correlated  having  Gaussian  power  spectrum 
centered  at  zero  frequency  with  a  standard  deviation  of  Op.  Thus,  the  desired 
power  spectrum  is 

4»ff(f)  =  —A -  exp(-f2/2oF2)  .  (G.l) 

\2n  OF 

Consider  a  continuous  linear  system  whose  impulse  response  is 

g(t)  -  exp(-t2/o2)  (G.2) 

where 


rT  =  l/2noF  . 


(G.3) 


The  transfer  function  of  this  non-causal  system  is 


G(f  ) 


— -I -  cxp(-f^/4ap)  . 

V^°F  F 


(G.4) 


Let  the  input  to  this  system,  denoted  n(t),  be  a  zero  mean,  unit  variance,  Gaussian 
noise  process  with  a  flat  spectrum. 


*nn(f> 


S  for|  f  |<  A_ 
2S 


otherwise 


(G.5) 
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Then,  the  output  process  c(t)  is  also  a  zero  mean  Gaussian  process  whose  power 
spectrum  is 


•  l«f>l2»  <» 


=  -5__  exp(-f2/2op) 
/.  _  -  ^  r 


for  f  <  4- 


otherwise 


(G.6) 


The  restriction  on  the  region  of  existence  can  be  removed  if  S  is  chosen  to  be 
small  enough  such  that 


. 


(G.7) 


From  its  power  spectrum,  it  is  obvious  that  if  the  input  noise  is  sampled  at 
intervals  of  S,  the  samples  are  independent  and  have  a  Gaussian  distribution  with 
zero  mean  and  unit  variance.  Thus,  a  discrete  equivalent  of  the  continuous  system 
may  be  written  as. 


c(t)  -  S  ]£n(kS)  g(t-kS) 
k 

If  c(t)  is  sampled  at  intervals  of  6, 
c(m5)  =  S  ^  n(kS)  g(mS  -kS) 


(G.8) 


=  S.  f(m)  . 


(G.9) 


The  power  spectrum  of  (f(m)}  is, 


«££(f)  -  *cc(0 

-  - - — ^  exp(-f2/2cjp) 

S.  4trap 

=  - -== - exp(-f2/2a2)  (G.10) 

2  Er»r  r 

where  S  *  "l  (G.ll) 

Note  that  the  definition  of  S  satisfies  the  inequality  given  in  (G.7)  and  that 
the  power  spectrum  of  {f(m)}  is  the  same  as  that  of  the  required  clutter 
components.  Thus,  the  clutter  sequence  can  be  generated  using 

f(m)  «  ^  n(kS)g(m6-  kS)  .  (G.12) 

k 

However,  it  must  be  noted  that  the  impulse  response  of  the  filter  has  infinite 
duration.  For  implementation  in  a  digital  simulation,  the  impulse  response  is 
truncated  to  1 1 1 <3S .  Since  the  input  noise  sampled  at  intervals  of  S,  six  samples 
of  noise  are  required  in  the  summation  in  (G.12).  Note  however,  that  several 
samples  of  clutter  can  be  generated  from  the  same  six  noise  samples  because  the 
output  sampling  interval  6  is  much  smaller  than  S.  In  fact,  I  samples  of  clutter 
can  be  computed  from  six  samples  of  input  noise  where  I  «=  S/6. 


APPENDIX  H 


PROGRAM  LISTINGS 


This  appendix  contains  the  listings  of  all  the  computer  programs  developed  for 
this  study.  The  programs  re  in  FORTRAN  and  are  compatible  for  execution  on  the 
CDC-6600  computer  system. 


i-UBROUTINE  ANIT  74/74  OPT=l 


SUBROUTINE  ANIT (START) 

^•initialization  for  noise  generators 

<;••••••••••••••••«••••«••*•««« . . 

5"  C0MM0N/RANDM/AR0*RN01 *ps 

AROsSTART 


RN01=START*0» OOOOOl 

RS*47436.0 

RETURN 

ENO 


FTN  4*2*74355  1 


U8R0UTINE  CLUTTN 


74/74  0PT=1 


FTN  4*2*74355 


SUBROUT IKE  CLUTTN ( X , SCR  * XM2 * SIGHAF , NCOEL > 
£***************»**»************«**•**«»***************•**•**********•** 
C**GENERATES  SAMPLES  OF  GROUNO  CLUTTER  AND  ADDS 
C*«THEM  TO  THE  COMPLEX  ARRAY  X 

c«*x  is  the  input/output  complex  array 

C*«SCR  IS  THE  SIGNAL  TO  CLUTTER  RATIO 
C**NOTEs=SIGNAL  power  is  assumed  to  be  unity 
C**XM2  IS  THE  OC  TO  AC  POWER  RATIO 
C**AC  CLUTTER  HAS  GAUSSIAN  SPECTRUM 
C*«WITH  ZERO  MEAN 

C**SIGMAF  IS  THE  STANDARD  DEVIATION  OF 
C**THE  CLUTTER  SPECTRUM 

C**NCDEL  IS  THE  RANGE  DELAY  CORRESPONDING 
C**TO  THE  CLUTTER  LOCATION 

£*#*«#******•»**** ************************** a**************************** 
COMMON/ADPAR/DELT *QAD*VMAX 

COMMON/PROPAR/NLOOK ,NSAMPP ,NFFT .NSCLK  *NC INT 
COMMON/XCOOE/NSTAGE ♦ INI T ( 1 0 ) »FEEDBK (10) *COOE ( 1024 ) 

COMPLEX  X<1) 

DIMENSION  XR (4) *XN ( 12) 

INTEGER  CODE 
C0NST=1 .0/ (SORT (SCR) ) 

NCYCLE=NLOOK/NSAMPP 
CALL  RAnDU(XR*2) 

PI=3. 141492654 
PSI=2.C*PI*XR(2) 

CPSI=COS(PSI) 

SPSI*SIN(PSI) 

IF(XM2.LT.O)  GO  TO  20 
G0=SQRT(XM2/(1.0*XM2) ) 

GA=SQRT(0.5/(1.0*XM2) ) 

GC*G0*CPSI 

GS*GO*SpSI 

CALL  RAnDG(XN*12*0. 0*1.0) 

PI2s2.0*PI 

SIGTAU=1.0/(PI2*SIGMAF) 

SIGSS=SIGTAU/DELT 
SIGI=SIGSS/1. 4142 13562 
CDEN=0.5/(SIG1*SIG1) 

JTIb1«253314137*SIGSS 

TI=FL0AT(ITI) 

00  10  IC=1 *NCYCLE 
DO  l l  10s 1 tNSAMPP 
J=(IC-1)*NSAMPPMD 
IC0DE=ID-NC0EL 

IF(ICODE.LE*0)ICODE=ICODE»NSAMPP 

CIsO.O 

CQ=0«0 

DO  12  1=1*6 

M*(6-I)*ITIO 

BsFLOAT (M) 

DUM«;B-3.0«TI 

EXPT=-DUM»DUM*CDEN 

HM=EXP(EXPT) 

CIsCI*XN(I)»HM 

12  CQ=CQ»XN(1*6)*HM  - 

CITsFLOAT (CODE ( ICOOE) > * (GC*GA*C1 ) 

H-3 


Subroutine  cluttn 


74/74  .  OPT=l 


60 


65 


70 


CQTsFLOAT (CODE ( ICOOE) > « <GS*GA*CQ> 
X(J)*X(J>*  CONST»CMPLX(CIT»COT) 

11  CONTINUE 
10  CONTINUE 
RETURN 
20  CONTINUE 

00  30  IC=1»NCYCLE 
00  31  I0=1*NSAMPP 
JMIC-1)«NSAMPPMD 
ICOOE=IO-NCOEL 

IFfICODE.LE.O)ICODE=lCODE*NSAMPP 
CITsFLOaT (CODE ( ICOOE) )«CPSI 
CQTsFLOAT (CODE (ICODE) >  *SPSI 
X(J)*X(J)*  CONST«CHPLX(ClT*CQT) 

31  CONTINUE 
30  CONTINUE 
RETURN 
ENO 


75 


SUBROUTINE  CORREL  74/74  OPT=l  FTN  4. 2*74355  1 


SUBROUTINE  CORREL <X*NlN»NOUT> 


C**PULSE  COMPRESSION  FILTER  FOR  THE 
*  COPROCESSORS  USING  PULSE  COMPRESSION 

5  C**X  IS  THE  COMPLEX  INPUT/OUTPUT  ARRAY 

r  C**NIN  IS  THE  NUMBER  OF  INPUT  SAMPLES 

C**NOUT  IS  THE  NUMBER  OF  OUTPUT  SAMPLES 

C************************** ************************************ *«•«««••• 

C0MM0N/PR0PAR/NL00K,NSAMPP,NFFT*NSCLK»NCINT 

10  COMMON/xCODE/NSTAGE* INIT (10) »FEEDBK(10>  »C0DE(1024> 

COMPLEX  X ( 1 > *XOUT 
INTEGER  CODE 
NOUT=NlN*NSAMPP 
NSMPP 1 =NSAMPP* 1 

15  NINP=NIN*1 

DO  10  II=NINP»NOUT 

I=NOUT-II*NINP 

XOUT=CMPLX(0.0«0.0) 

JFIN=N0UI-I 

20  DO  11  J=1 ♦ JFIN 

IND=J* I-NSAMPP 

IF (CODE (J) «GT«0)G0  TO  12 

XOUT=XOuT-X ( IND) 

GO  TO  11 

25  12  XOUT-XOuT *X ( IND) 

11  CONTINUE 

X ( I ) =XOuT 
10  CONTINUE 

DO  20  I I=NSMPP1 » NIN 

30  I=NIN-I I ♦NSMPP 1 

XOUT=CMPLX(O.OtO.O) 

DO  21  J=1»NSAMPP 

IND=J*I-NSAMPP 

IF (CODE (J)*GT«0)G0  TO  22 

35  X0UT=X0UT-X(IND) 

GO  TO  21 

22  XOUT=XOuT*X (IND) 

21  CONTINUE 

X(I)=X0UT 

40  20  CONTINUE 

DO  30  1 1  =  1 *NSAMPP 
I=NSAKPp-I 1*1 
XOUT=CMPLX(0.0»0.0) 

JST=NSAMPP-l*i 

45  DO  31  J=JST»NSAMPP 

INO= J* I -NSAMPP 

IF (CODE ( J) *GT«0)G0  TO  32 

X0UT=X0UT-X (IND) 

GO  TO  31 

50  32  X0UT=X0UT*X<IND> 

31  continue  y 

X  ( I )  =X0ljT  | 

30  CONTINUE 

RETURN 

55  ENO  ij 


f  SUBROUTINE  OECOOE 


74/74  OPT=l 
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SUBROUTINE  OECOOE (X*N»IR) 

<;««•••»««••••»»**»»•»••» . 

c*« IMPLEMENTS  THE  DECODER 

C*«X  IS  THE  COMPLEX  INPUT/OUTPUT  ARRAY 

C**N  IS  THE  NUMBER  OF  SAMPLES 

C*«IR  IS  THE  REQUIRED  RANGE  BIN 

£******»«*«*»****»***»#**«#*********«*•*****»***»*«*«***«******«***»**** 
COMMON/PROPAR/NLOOK ,NSAMPP .NFFT ,NSCLK  tNC INT 
COMMON/xCODE/NSTAGE ♦ INI T ( 10 ) .FEEOBK ( I 0 ) »COOE (1024) 

COMPLEX  X(l) 

INTEGER  CODE 
NCYCLE=N/NSAMPP 
DO  10  ICS1 .NCYCLE 
00  10  I0=1*NSAMPP 
J=(IC-1)*NSAMPP*ID 

icooe=id-ir 

IF(ICODE.LE.O)ICODE=ICODE*NSAMPP 
IF (CODE ( ICODE) .LT.O) X ( J) =-X ( J) 

10  CONTINUE 
RETURN 
END 
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55 


SUBROUTINE  DFT (A.N* ISN.NRlN) 

C »#*****#**#**#*****#«#»**«**«**«*«*«*»**««*•*•*»#*«*****#* »## ********** 

C**COMPUTES  The  DISCRETE  FOURIER  TRANSFORM  OF 

C*«THE  COMPLEX  ARRAY  X 

C*«X  IS  THE  COMPLEX  INPUT/OUTPUT  ARRAY 

C*«N  IS  THE  NUMBER  OF  SAMPLES  IN  X 

C**ISN  =  -1  FOR  DIRECT  DFT 

C**ISN  =  *1  FOR  INVERSE  DFT 

C**NB1N  IS  NUMBER  OF  OUTPUT  SAMPLES  DESIRED 

C**IF  N  =  2*»M  ♦  FFT  IS  USED 

C**0THERWISE  DFT  is  used 

£«««««  *#<******■««•&**«**•#**«*«***  ****#«*«**«**«***«*••**««*«•***»****«««*« 
COMPLEX  A  < 1 )  ♦T1.T2.TEMP,X(256) 

PI 2=6 .283 18530 7 17959 

MASK=1 

IGAM=0 

1  IGAM=IGAM*1 
M=2**IGAM 
IF(N-M)2,4»1 
4  CONTINUE 
N1=N-1 

DO  20  I I=2.N1 
1=11-1 
IFLIP=0 

DO  10  J=1 . IGAM 
LSN=2**(J-1) 

10  IFLIP=2*IFLIP* (MASK. AND. ( I/LSN) ) 

IF (I.LE.IFLIP)GO  TO  20 
11=1*1 
I2=IFLIP* 1 
TEMP=A  < 1 2 ) 

A(I2)=A(I1) 

A(I1)=TEMP 
20  CONTINUE 

DO  30  1=1. IGAM 
NEL=2**I 
NEL2=NEL/2 
nset=n/nel 

ANG=PI2/NEL 

SI=SIN(ANG) 

CI=COS<ANG) 

DO  30  J=1.NSET 
INCR=<J-1)*NEL 
S0=0.0 
C0=1.0 

DO  30  1 1  =  1 .NEL2 
J1=II*INCR 
J2=J1*NEL2 
T1=A(J1) 

T2=A(J2)»CMPLX(C0.ISN«*S0) 

A<Jl)=Tl*T2 
A(J2)=T1-T2 
SN=S0*CI*C0«SI 
CS=C0*CI-S0»SI 
C0=CS 
30  SO=SN 

GO  TO  100 
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50 


55 


70 


75 


2  CONTINUE 

ANG=ISN*PI2/FL0AT (N) 

DO  40  I=l*NBIN 
X(I)=CMPLX(0. 0*0.0) 

00  50  J=1*N 

ARG=ANG*FLOAT ( ( 1-1 ) * ( J-l ) ) 
T1=CMPI_X(0.0*ARG) 

T2=CEXP (T  I ) 
X(I)=X(I>4T2*A(J) 

SO  CONTINUE 
40  CONTINUE 

DO  60  1=1 *N 
60  A  ( I )  =X  ( I ) 

100  CONTINUE 

IF(ISN.GT.O)RETURN 
DO  110  1=1 *N 
A(I)=A(I) /FLOAT (N) 

110  CONTINUE 
RETURN 
END 
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SUBROUTINE  FFTOEC (X«NSTART *NSKlP*NBIN»NOUT) 
c*******«*«***«*****»*****«*******************»******************«****** 
C»»FFT  PROCESSING  FOR  THE  OECODER  APPROACH 
C**X  IS  THE  COMPLEX  INPUT/OUTPUT  ARRAY 
C*«NSTART  IS  SAMPLE  WHERE  INTEGRATIN  BEGINS 
C**NSKIP  SAMPLING  RATE  REDUCTION  FACTOR 
C*«NBIN  IS  THE  NUMBER  OF  OUTPUT  PINS  DESIRED 
C*«NOUT  =  NBIN  IN  THIS  SIMULATION 

•*»**•*****• *««***********«**#****»***««***«***#**«***<»*************«. 
COMMON/PROPAR/NLOOK ,NSAMPP  *  NFFT ,NSCLK  »NC INT 
COMPLEX  X ( 1 ) tXFFT (2561 
NS=  1 

ISTIN=NSTART 

ISTOUT=0 

NOUT=NBlN*NCINT 

DO  10  INT=1,NCINT 

DO  11  IfFT=1*NFFT 

ISAMP=(IFFT-1)«NSKIP*ISTIN 

XFFT(IFFT)=X<ISAMP) 

11  CONTINUE 

CALL  WEIGHT (XFFT »NFFT ) 

CALL  DFt<XFFT*NFFT,-1»N8IN) 

DO  12  IfFT=1 *N8IN 
X(IFFT*IST0UT)=XFFT(IFFT> 

12  CONTINUE 
ISTOUT=lSTOUT ♦NBIN 
ISTIN=ISTIN4NFFT«NSKIP 

10  CONTINUE 
RETURN 
END 
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SUBROUTINE  FFTPC (X*NSTART+NBIN> 

##•*****«**##* #**#**•*#•**•*»«•*******« 

C**FFT  PROCESSING  FOR  PULSE  COMPRESSION 
^•CONFIGURATIONS 
C**X  IS  THE  COMPLES  INPUT  ARRAY 
C**NSTART  SAMPLE  AT  WHICH  INTEGRATION  STARTS 
C**NBIN  IS  THE  NUMBER  OF  OUTPUT  BINS  DESIRED 
(««»»«••*«•>««*»«+««»»•««««*««««««««•«<»«««•««•«»«««««*•#«««••»«««••««««{ 
C0MM0N/PR0PAR/NL00K,NSAMPP*NFFT*NSCLK»NCINT 
COMPLEX  X(l) ♦XFFT(256) 

NCOH=NFFT«NSAMPP 

NFIN=NSTART+ (NCOH«NCINT)-I 

DO  10  int=nstart»nfin*ncoh 

DO  11  ISAMP=1*NSAMPP 

DO  12  IFFT=1«NFFT 

J=(INT-1 >+<IFFT-l)«NSAMPP+ISAMP 

XFFT(IFFT)=X(J) 

12  CONTINUE 

CALL  WEioHT  <XFFT *NFFT) 

CALL  DFt(XFFT»NFFT*-1*NBIN) 

DO  13  IFFT=1*NFFT 

J= ( INT-1 ) ♦ ( IFFT-1 ) *NSAMPP+ ISAMP 

X(J)=XFFT(IFFT) 

13  CONTINUE 
11  CONTINUE 
10  CONTINUE 

RETURN 

END 
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SUBROUTINE  FILCAS(X,NSTART*NFIN> 

. . . . . . . 

C*MMPLEMENTS  RECURSIVE  OISCRETE  FILTERS 
C**AS  CASCADE  OF  SECOND  AND  FIRST  ORDER  SECTIONS 
C**CANONICAL  REALIZATION  IS  USED 
C**X  IS  THE  COMPLEX  INPUT/OUTPUT  ARRAY 
C*»NSTART  IS  SAMPLE  WHERE  FILTERING  BEGINS 
C**NF IN  IS  SAMPLE  WHERE  FILTERING  ENDS 

<;•«••••««•«««»«»«««««•*««» . <»••«••••••••••••«•• . . 

COMMON/FLTT/NSEC*NFIRST  » NTYPE (10) *CSEC<2*10) »CFIRST. CONST 

COMMON/PROPAR/NLOOK,NSAMPP,NFFT,NSCLK*NCINT 

LOGICAL  NTYPE 

COMPLEX  X  < 1 > ,Y<3) ,DUM 

DO  10  IsEC=l*NSEC 

DO  11  1=1*3 

11  Y(I)=(0. 0*0*0) 

DO  12  I=NSTART*NFIN 
Y(l)sX(I) 

DO  13  I F - 1 »  2 

13  Y(l)=Y<l>-CSEC<IF*ISEC>*Y(IF*l> 

DUM=Y ( 1 ) *Y (3) 

IF (NTYPE ( ISEC) ) GO  TO  15 
DUM^OUM-Y (2)-Y (2) 

GO  TO  16 

15  DUM=DUM.Y (2) *Y (2) 

16  CONTINUE 
X(I>=DUm 
Y (3) =Y (?) 

Y (2  >  =Y  ( 1 ) 

12  CONTINUE 
10  CONTINUE 

IF(NFIRST.EQ.0)G0  TO  25 
Y(l>=(0. 0,0.0) 

Y(2)=(0. 0,0.0) 

DO  22  I=NSTART ,NFIN 
Y  ( 1 ) =X ( I ) 

Y ( 1 ) =Y ( 1 ) -CFIRST»Y (2) 

IF(NTYPE(NSEC*1)  )G0  TO  23 
DUM=Y(1)-Y(2) 

GO  TO  24 

23  DUM=Y ( 1 ) ♦ Y (2) 

24  CONTINUE 
X(I)=DUM 
Y (2) =Y ( 1 ) 

22  CONTINUE 

25  CONTINUE 

DO  30  IsNSTART *NFIN 
30  X(I>=X(I)*CONST 
RETURN 
END 
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SUBROUTINE  INTPC (XtNSTART) 

. . . . 

C**NON  COHERENT  INTEGRATION  FOR  THE 
C*»PULSE  COMPRESSION  APPROACH 

*************************** *•*•*****••*•*•*•** 

COMMON/pROPAR/NLOOK,NSAMPP*NFFT.NSCLK*NCINT 

COMPLEX  X ( 1 ) *X INT 

NCOH=NFFT  *nsampp 

00  10  1=1 »NCOH 

XINT=CMPLX(0. 0*0.0) 

00  11  int=i.ncint 

J= (NSTART-1 ) *1 ♦ ( INT-1 ) *NCOH 
XINT=XInT4X(J) 

11  CONTINUE 
X<I>=XInT 
10  CONTINUE 
RETURN 
END 
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SUBROUTINE  MTI <X«NIN*NPER*NPULSE*NOUT> 

. . * . . . . 

COMPLEMENTS  THE  DELAY  LINE  CANCELLER 
C»*X  IS  THE  COMPLEX  INPUT/OUTPUT  ARRAY 
C«*NIN  IS  THE  NUMBER  OF  INPUT  SAMPLES 
C*«NPULSE  IS  NUMBER  OF  PULSE  CANCELLED 
C*«NPER  IS  NUMBER  OF  CODE  PERIODS  IN  DELAY 
C**NOUT  IS  NUMBER  OF  OUTPUT  SAMPLES 


COMMON/pROPAR/NLOOK  *NSAMPP»NFFT  »NSCLK*NCINT 
COMPLEX  X ST ORE ( 150 ) « YSTORE ( 1 50) *X  ( 1 ) 

MlNsNIN 

noelay=nsampp*nper 
DO  15  IFILT=1*NPULSE 
MOUT=MlN*NDELAY 
DO  10  1  =  1 »NDELAY 
XSTORE ( I ) =X  ( I ) 

10  CONTINUE 
NDEL1=N0ELAY*1 

DO  11  I=NDEL1*MIN»NDELAY 
DO  12  J=1 »NDELAY 
YST0RE(J)=X(I*J-1) 

12  CONTINUE 

DO  13  J=1 » NDELAY 
X ( I ♦ J-l ) = YSTORE (J) -XSTORE ( J) 

13  CONTINUE 

DO  14  Jsl »NDELAY 
XSTORE ( J) = YSTORE ( J) 

14  CONTINUE 

11  CONTINUE 
MINP=MIN*1 

00  16  J=MINP,MOUT 
X(J)=-YST0RE(J-MIN) 

16  CONTINUE 
MINsMOUT 

15  CONTINUE 
nout=mout 

RETURN 

END 


I 
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SUBROUTINE  NOISE (X) 

C**ADDS  NOISE  TO  THE  INPUT  SIGNAL 
C*»X  IS  THE  COMPLEX  INPUT/OUTPUT  ARRAY 

. . . . * . * . . 

COMHON/xNOISE/XMEAN*STOEV 

COMMON/pROPAR/NLOOK ♦NSAMPR»NFFT*NSCLK  «NC INT 
DIMENSION  XU (24) 

COMPLEX  X ( 1 ) *GAUSS 
DO  10  1=1 *NLOOK 
CALL  RAnDU (XU»24) 

GAUSS1=0*0 
GAUSS2*0.0 
DO  11  J=1 t 12 
GAUSS1=GAUSS1*XU(J) 

11  GAUSS2*GAUSS2*XU(J*12) 

GAUSS1=(GAUSS1-6.0)»STDEV»XMEAN 
GAUSS2=(GAUSS2-6.0)«STDEV^XMEAN 
GAUSS=Cfc.JLX  (GAUSS1  «GAUSS2> 

10  X(I)=X(I> ♦GAUSS 

RETURN 
END 
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SUBROUTINE  PLOT  IP (X*MX,MS*IPX*MPtMSI7E*FRM»KPERI0) 


••GENERATES  A  LINE  PRINTER  PLOT 


DIMENSION  X(l)  tXL(UO)  «W(6) 

DATA  BLANK *ASTR I X «HOR  t VERT/ 1H  ♦ IH*» lH-» IHI/ 
DATA  W(1)/10H(1X#I5.2X,/ 

DATA  W ( 3) /3HA1 */ 

DATA  W (4) /9HlX*EI2»6) / 

PERIODIC  PLOT 

X*ARRAY  TO  BE  PLOTTED 

MX*STARTING  ADDRESS 

HSKIPX=DISTANCE  BETWEEN  POINTS 

MP  *  NUMBER  OF  POINTS  plotted 

MSIZE  =  WIDTH  OF  THE  PLOT  (LESS  THAN  1 10 > 

FRM  *  "MSIZE" 

KPERIO  =  PERIOD 
WRITE (6*1) 
l  FORMAT ( 1H1 ) 

IF(MSIZE.GT.llO)  MSIZE=ilO 

IF(MSIZE.GE,1I0)FRM=3H110 

W  (2) =FRm 

MMX«MX»mP-1 

MA=MX 

IF(MA.LF.O)  MAaMA^KPERIO 
XMIN=X(MA) 

XMAX=X(mA> 

DO  5  JJsMX*MMX*MSKIPX 
JsJJ 

IF(J.LE.O)  J=J*KPERIO 
IF(J.GT. KPERIO)  J=J-KPERI0 
IF(X<J> .GT.XMAX)  XMAX=X(J) 

IF(X(J) .LT.XMIN)  XMIN=X(J) 

5  CONTINUE 

IF( (XMAX-XMIN) .LE.l.E-7>XMlN=XMIN-l.E-5 
KS=  1 

IF(XMIN.GT.O.O.OR.XMAX.LT.O.O)  ks=? 
DELXs(XMAX-XMIN) /FLOAT  (MSlZE-1) 

GO  TO  <7«8) *KS 
7  NAXIS=1.0-XMIN/DELX 
8  CONTINUE 

DO  13  JA=MXtMMX.MSKIPX 
JXsJA 

IF(JX.LE.O)  jx*jx»kperio 
IF(JX.GT. KPERIO)  jx=jx-kperio 
DO  14  J=l* MSIZE 
14  XL(J)=BLANK 

IF(KS.EO.l)  XL(NAXIS)*VERT 
IX*1.0^(X(JX)-XMIN)/DELX 
IFMX.GT. MSIZE)  I  X=MS  I ZE 
IF(IX.LT.l)  I X» 1 
GO  TO  (2l*22)tKS 
21  IF(IX.LF.NAXIS)  GO  TO  23 
NAXISP=NAXIS*1 
DO  24  I =NAXISP* IX 
24  XL ( I ) ®AsTRIX 
23  IFUX.NF.NAXIS)  GO  TO  25 
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XL(IX)=aSTRIX 

60  TO  26 

6C 

25 

NAXISM=nAXIS-1 

00  27  I=IX*NAXISM 

27 

XL(I)=AsTRIX 

GO  TO  26 

22 

!F(XHIN.GT.O.)  GO  TO  30 

65 

DO  31  IslXtMSIZE 

31 

XL ( I ) sAsTRIX 

GO  TO  26 

30 

DO  28  I =1 • IX 

28 

XL(I)=AsTRIX 

7C 

26 

continue 

WRITE (6 tW) JAt (XL (I) » 1=1 .MSIZE > »X < JX> 

13 

CONTINUE 

RETURN 

END 

•2*74355 
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SUBROUTINE  PNCODE 

£*•»« #****#*#*»**#**##*•»****#»•«#»*** •  *«■»**#**#****#****##***«***#**#*< 

C**GENERATION  OF  THE  PSEUOO  RANDOM  SEQUENCE 
£»*******»**•*»********** #»*•**#«»***#***#*###*»#*»**##***«*»#*#*****#** 
C0MM0N/PR0PAR/NL00K,NSAMPP,NFFT*NSCLK»NCINT 
COMMON/xCODE/NSTAGE « INI  T (10) »FEED8K (10) »CODE  <1024) 

DIMENSION  RGSTR (10) 

INTEGER  F£ED8K*C0DE, RGSTR 
DO  10  Isl *NSTAGE 
10  RGSTR ( I )  =  INIT  < I ) 

DO  1  IC0DE=1*NSAMPP,NSCLK 
KBACK=0 

DO  2  I=1*NSTAGE 

KBACK=KBACK*FEEDBK ( I ) *RGSTR ( I ) 

2  CONTINUE 

KBACKH=KBACK/2 

icodep=icode*nsclk-i 

DO  4  JC0DE=IC0DEtIC00EP 
CODE(JCODE)=RGSTR(NSTAGE) 

4  CONTINUE 

DO  3  J«?»NSTAGE 

I=NSTAGe*2-J 

RGSTR ( I ) =RGSTR  1 1-1 ) 

3  CONTINUE 

RGSTR  < 1 ) =KBACK-2*KBACKH 
1  CONTINUE 
RETURN 
END 


I 
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J  SUBROUTINE  QUANT (X) 

£*****»*« «•«»•«««••««•»•«««»««••••«»•«»«•«•••«•««««••«««*«•*••««« 0 «««»•« 
I  C*»QUANTIZATION  =*  a/o  conversion 

j  c •**«•**** « ***»*«*********««**«*#**»*«**«•***»# «*«««•««« «•••«««•« ««««•«« 

5  COMMON/ADPAR/D£LT*OAD»VMAX 

.  COMMON/PROPAR/NLOOK.NSAMPP,NFFT.NSCLK»NCINT 

COMPLEX  X(l) 

*  00  10  1=1 »NLOOK 

X(I»=X<I)/VMAX 

10  XA=REAL(X(I>) 

XB=AIMAG(X(I) ) 

AX=A8S (XA) 

BX=ABS(XB) 

I AX=AX/QAO 

15  !BX=BX/QAD 

AY=QAD*(FL0AT(IAX)*0.5) 

BY=QAD* (FLOAT (IBX) *0.5) 

XI=SIGN(AY»XA) 

XQ=SIGN (BYtXB) 

20  X  ( I )  =CMPLX  (XI  *  XQ) 

10  CONTINUE 
RETURN 
END 


I 
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SUBROUTINE  RANOG (X«N»XM£AN*STDEV) 

•»***#*#**** » «*«»«** « •««»«* »***«#******#********« ************* 

C**GENERATES  GAUSSIAN  RANDOM  NUMBERS 
C*«N  IS  THE  NUMBER  OF  NOISE  SAMPLES 
C**X  IS  THE  OUTPUT  ARRAY 

C**XMEAN  IS  THE  MEAN  OF  THE  NOISE  SEQUENCE 
C**STDEV  IS  THE  STANDARD  DEVIATION 
C**OUTPUT  IS  AN  INDEPENDENT  SEQUENCE 

C*«************ *********** s********************************************* 

DIMENSION  X ( 1 ) t XU ( 1 2 ) 

DO  10  1  =  1  *N 
CALL  RANDU<XU»12) 

GAUSS=0.0 
DO  11  J=l,12 
11  GAUSS=GAUSS*XU(J) 

GAUSS=(GAUSS-6.0)»STDEV*XMEAN 
10  X(I)=GAuSS 
RETURN 
END 


1 
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'  SUBROUTINE  RANDU(X.N) 

C## »*»**« #«*«*»»»**» ««««««««« **»**»#***#*«*#*#****«***•***##*  *##*****«»** 

j  C**GENERATES  UNIFORMLY  DISTRIBUTED 

I  C**RANDOM  NUMBERS 

5  C*«N  IS  THE  NUMBER  OF  NOISE  SAMPLES 

[  C*«X  IS  THE  OUTPUT  ARRAY 

I  C««OUTPUT  IS  AN  INDEPENDENT  SEQUENCE 

COMMON/RANDM/ARD.RND1 »RS 
10  DIMENSION  x<l) 

DO  3  1=1 »N 
ARN=ARD#«2*RS**2 
K=ARN/1 OOOOOOOOO* 

ARD5(ARN-FLOAT(K)*l OOOOOOOOO. )/l 000. 

15  IF ( ARD) ?* 1 *2 

•  1  ARD=1.0 

2  RS=RS* 1 
RND=ARD*0. 000001 

3  X(I)=RND 

20  RETURN 

END 
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SUBROUTINE  RMS (X.NSTART ) 

£««•«• **«•»**» «»««««« *«***» »» » ******************* «««««»« 

C«*COMPUTES  THE  RMS  value  of  each  sample  of 

C*«THE  COMPLEX  ARRAY  X 

C««RESULT  IN  REAL  PART  OF  ARRAY  X 

COMPLEX  X(l) 

COMMON/pROPAR/NLOOK.NSAMPPfNFFT .NSCLK  t NCI NT 

ncoh=nfft*nsampp 

NFIN=NSTART* (NCOH*NCINT)-l 
DO  10  I =NSTART *NF IN 
SSQ=CABS(X(I) ) 

X ( I ) =CMPLX (SSQtO.O) 

10  CONTINUE 
RETURN 
END 
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SUBROUTINE  SIGNAL (X) 

C**GFN£RATES  SAMPLES  OF  SYNTHETIC  VIDEO 
C**INCLUDES  TARGET  RETURN^CLUTTER  AND  NOISE 
C**CLUTTER  INCLUDES  IN  RANGE  CLUTTER 
C*«OUT  OF  RANGE  DISTRIBUTED  CLUTTER 
C**FIX£D  CLUTTERf  ANTENNA  SPILLOVER 
C**AND  OTHER  ISOLATED  CLUTTER 

COMPLEX  X ( 1 > 

COMMON/S IMPAR/NMONTE « INOPT (7) * IPRINT (10) 

COMMON/pROPAR/NLOOK,NSAMPP,NFFT,NSCLK,NCINT 

COMMON/SIGPAR/FDOP.NDELTA 

COMMON/CL TP AR/SCR I *SCRD « SCRF  »SPIL  *SIGMAF *XM2»NCDELA*XM20 t SIGMAO 
*  fNCDELOtSCRO 
DO  10  1=1 »NLOOK 

10  X<I)=CMPLX(0.0»0.0) 

IF (INOPT  < 1 ) »EQ. 1 ) CALL  TARGET(X) 

IF(IN0PT(2)  .EQ.DCALL  NOISE  (X) 

IF(INOPT(3)*EQ.O)GO  TO  11 

CALL  CLUTT(X*SCRI*XM2»SIGNAF*NDELTA) 

11  IF(IN0PT(4) ,EO.O)GO  TO  1? 

DO  13  ICELL=1»NSAMPP,NSCLK 
IRANG= ( ICELL- 1 ) /NSCLK 
DENOM=1.0*0.15*FLOAT(IRANG> 

WEIGHT=(1.0/DEN0M)««3 
SCRDR=SCRD«3. 4/WEIGHT 

CALL  CLUTT(X»SCRDR,XM2,SIGMAF,ICELL) 

13  CONTINUE 

12  IF ( INOPT (5) *EQ« 0 ) GO  TO  14 

CALL  CLUTT ( X  4  SCRF ♦ - 1 . 0 , 0 . 0 , NCDELA ) 

14  IF (INOPT (6) *EQ»0)G0  TO  15 

CALL  CLUTT (X,SCR0»XM20,SIGMA0»NCDEL0) 

15  IF ( INOPT (7)*EQ»0)G0  TO  16 
CALL  CLUTT (X»SPIL*-1.0«0.0t0) 

16  CONTINUE 
RETURN 
END 
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SUBROUTINE  SPCTRM (X.N.NPLOT. IOPT) 

C««COMPUTES  ANO  PLOTS  THE  SPECTRUM  OF  X 
C«*X  IS  THE  COMPLEX  INPUT  ARRAY 
C**N  IS  THE  NUMBER  OF  SAMPLES  IN  X 
C**IOPT  =  l  FOR  MAGNITUDE  PLOT  ONLY 
C**IOPT  =  2  FOR  PHASE  PLOT  ONLY 
C**IOPT  =  3  FOR  BOTH  PLOTS 

C ********** »»»«»«««« tt »««»»« » «««•*•«»»«« »««*»»«« »»•!»«»«•««»««« ««««««« $4 «« 

COMPLEX  X  (1 ) *Y (1024) 

DIMENSION  PLOT (1024) 

DATA  FRM/3H100/ 

DO  10  1=1. N 
10  Y(I)-X(I) 

CALL  DFT(Y.N.-l.N) 

Pl=3. 141592654 

Pl2=PI/?*0 

NP2=NPL0T/2 

NP2P=NP2*l 

NP2M=NP2-1 

IF(I0PT.EQ.2)G0  TO  15 
DO  20  1=1 *NP2P 

20  PL0T(I*NP2M)=CABS(Y(D) 

DO  21  I=1.NP2M 

21  PLOT ( I ) =CABS (Y (N*I-NP2M) ) 

CALL  PL0T1P(PL0T.1.1.NPL0T.100.FRM.NPL0T> 

IF { IOPT. EQ. 1 ) RETURN 
15  CONTINUE 

DO  30  1=1 *NP2P 
YR=REAL (Y  ( I ) ) 

YI=AIMAG<Y(I) ) 

IF(ABS<YR) *LT.1.0E-09)G0  TO  31 
PLOT ( I *NP2M) =ATAN2 (YI »YR) 

GO  TO  30 

31  PL0T(I*nP2M)=SIGN(PI2»YI) 

30  CONTINUE 

DO  40  1=1 »NP2M 
YRsREAL (Y (N.I-NP2M) ) 

YI=AIMAG(Y<N*I-NP2M) ) 

IF(ABS(YR) .LT.1.0E-091G0  TO  41 
PLOT ( I ) =ATAN2  < YI »YR) 

GO  TO  40 

41  PLOT (I)=SIGN(PI2.YI) 

40  CONTINUE 

CALL  PLOT  IP (PLOT .1.1 .NPLOT  « 100.FRM.NPL0T) 

RETURN 

END 
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SUBROUTINE  START 

C*«« . **•* . . . 

(^•INITIALIZATION  SUBROUTINE 

C«*ALL  THE  PROGRAM  INPUTS  ARE  READ  IN  NAMELIST 
C**ARD.RND1»RS  ARE  INITIALIZATION  PARAMETERS  FOR 
C*«THE  NOISE  GENERATORS 

C**NBIT  IS  THE  WORD  LENGTH  OF  A/O  CONVERTER 
C*»FSAMP  IS  THE  SAMPLING  FREQUENCY 
C«*NSCLK  s  1  IN  THIS  SIMULATION 

C**VMAX  IS  HALF  THE  DYNAMIC  RANGE  OF  A/D  CONVERTER 

C*«SNRDB  IS  SlGNAL*TO»NOISE  RATIO  IN  DB 

C**NSEC  IS  NO.  OF  SECOND  ORDER  SECTIONS  IN  THE 

C««RECURSIVE  NOTCH/LP  FILTER 

C**NF IRST  IS  1  IF  ORDER  OF  FILTER  IS  ODD 

C*«OTHERWISE  nFIRST  IS  ZERO 

C**NTYPE  TYPE  OF  EACH  SECOND  OROER  AND  FIRST 

C**ORDER  SECTION  LOGICAL  ARRAY 

C**.TRUE.  FOR  LP  AND  .FALSE.  FOR  HP  FILTER 

C**NOTE==ONLY  BUTYERWORTH  AND  CHEBYSHEV  FILTERS 

C**CAN  BE  IMPLEMENTED  IN  THIS  SIMULATION 

C*«FDOP  IS  THE  DOPPLER  FREQUENCY  OF  TARGET 

C**NDELTA  IS  THE  TARGET  RANGE  BIN  NUMBER 

C*«NST AGE  IS  NO.  OF  STAGES  IN  THE  PN  SEQUENCE  GENERATOR 

C«*INIT  IS  THE  INITIAL  CONTENTS  OF  SHIFT  REGISTER 

C**FEEDBK  IS  THE  FEEDBACK  CONNECTION  FOR  THE 

C**SHIET  REGISTER  GENERATOR 

C**NFFT  IS  THE  NUMBER  OF  COHERENT  INTEGRATION  SAMPLES 

C*«NCINT  *  1  IN  THIS  SIMULATION 

C**NMONTE  IS  NOT  USED  IN  THIS  SIMULATION 

C««INOPT  IS  THE  INPUT  OPTION  ARRAY 

C**INOPT ( 1 ) =1  =  TARGET  IS  PRESENT 

C*«IN0PT(2)=1  =  NOISE  IS  PRESENT 

C**INOPT (3) =1  =  IN  RANGE  CLUTTER  PRESENT 

C*MNOPT  (4)  =1  =DISTRIBUTED  CLUTTER  IS  PRESENT 

C«*IN0PT(5)=1  =  FIXEO  CLUTTER  1$  PRESENT 

C*«INOPT(6)=l  =  OTHER  CLUTTER  IS  PRESENT 

C*MNOPT(7)=I  =ANTENNA  SPILLOVER  IS  PRESENT 

C**IPRINT  IS  NOT  USED  IN  THIS  SIMULATION 

c*«...m.scr  STANDS  FOR  SIGNAL  TO  CLUTTER  RATIO""" 

C**SCRD6I  IS  SCR  FOR  IN  RANGE  CLUTTER  IN  DB 

C*«SCRDBO  IS  SCR  FOR  DISTRIBUTED  CLUTTER  IN  OB 

C«*SCRDBF  IS  SCR  FOR  FIXED  CLUTTER  IN  DB 

C*«SCRDBO  IS  SCR  FOR  OTHER  CLUTTER  IN  DB 

C**SPILDB  IS  ANTENNA  SPILLOVER  IN  DB 
C««NCDELA  IS  RANGE  BIN  OF  FIXED  CLUTTER 
C*«SIGMAF  IS  THE  CLUTTER  SPREAD  FOR  IN  RANGE 
C**AND  DISTRIBUTED  CLUTTERS 

C«*XM2  IS  THE  CLUTTER  DC  TO  AC  ROWER  RATIO  FOR 

C»»IN  RANGE  AND  DISTRIBUTED  CLUTTER 

C»*NCDELO  IS  RANGE  BIN  OF  OTHER  CLUTTER 

C**SIGMAO  IS  SPREAD  OF  OTHER  CLUTTER 

C**XM20  IS  DC  TO  AC  POWER  RATIO  OF  OTHER  CLUTTER 

. . * . . **• . . 

COMMON/aDPAR/DELT»GAO»VMAX 
COMMON/PANDM/ARD»RNDI »RS 
COMMON/XNOISE/XMEAN.STDEV 

COMMON/FLTT /NSEC «NF IRST »NTYPE ( IO>  *CSEC (2*10) *CFIRST tCONST 


1 


»« 


(a 
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COMMON/ xCODE/NST AGE, INIT { 10> *FEEOBK<lO) *COOE<1024> 

,  COMMON/pROPAR/NLOOK  *NSAMPP  *NFFT ,NSCLK*NCINT 

60  common/sigpar/foop*noelta 

1  COMMON/SIMPAR/NMONTE*INOPT<7) *IPRINT  <IO) 

COMMON/CLTPAR/SCRI*SCRD«SCRF*SPIL*SIGMAF*XM2*NCOELA.XM20»SIGMAO 
*  *NCOELO*SCRO 
INTEGER  FEEDBK,CODE 
65  LOGICAL  NTYPE 

NAMELIST/RNDGEN/ARD,RND1 *RS 

NAMELIST/AD/NBIT*FSAMP,NSCLK*VMAX 

NAMELIST/NOIS/SNROR 

NAMELIST/FLT/NSEC,NFIRST*NTYPE*CSEC*CFIRST*CONST 

70  namelist/sig/fdop,ndelta 

NAMELIST/PNC/NSTAGE*1NIT*FEEDBK 

NAMELlST/PRO/NFFT*NCINT 

namelist/sim/nmonte*inopt*iprint 

NAMELIST/CLT/SCROBI*5CROBO,SCROBF*SCROBO*SPILDR*NCDELA,SIGMAF. 
75  *  XM2«NCDELO*SIGMAO*XM20 

NAMEUST/PRINT/N8I T  ,FSAMP  *NSCLK  ,  VMAX ,  SNRD8  »FDOP  ,NDELTA ,NST AGE » 

1  INIT*FeEOBK*NFFT*NCINT,NMONTE*INOPT,IPRINT*SCRDBI,SCRDBD* 

2  SCR0BF,SCR080,SPIL0B*NCDELA*SIGMAF,XM2,NCDEL0*SIGMA0*XM20 
3, NSEC, NFIRST*NTYPE.CSEC, CFIRST, CONST, ARD.RND1.RS 

80  READ <5* AC) 

REA0(5,RN0GEN) 

READ (5*nOIS) 

READ (5»FLT) 

READ (5,SlG> 

85  READ (5, PNC) 

READ (5*PR0) 

READ (5  *SlM) 

READ (5*CLT) 

WRITE(6, PRINT) 

90  delt=i.o/fsamp 

OAOsl.O/(2**(NBIT-l) ) 

SCRI=10««<0.1*SCRDBJ) 

SCRD=I0*«<0.1*SCRDBD) 

SCRF=10**(0.1*SCRDBF) 

95  SCRO=10««<O.I*SCRDBO) 

SPIL=10«»<0,1*SPILDB) 

SNR=10*»(0.1*SNRD8) 

VARsl ,0/SNR 
STDEV-SQRT (VAR/2. O) 

00  XMEAN=0.0 

LENGTHS (2*»NSTAGE)-1 
NSAMPPsLENGTH«NSCLK 
NCYCLE=NFFT*NCINT»5 
NLOOK*NCYCLE*NSAMPP 
05  CALL  PNCODE 

DO  10  I=1*NSAMPP 
10  CODE ( I ) *2*COOE ( I ) -1 
RETURN 
ENO 


] 
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SUBROUTINE  TARGET (X) 

. . . . * . * . 

C**AODS  SAMPLES  OF  TARGET  RETURN  TO  INPUT  ARRAY 
C*«x  IS  THE  COMPLEX  INPUT/OUTPUT  ARRAY 
C*«TARGET  IS  ASSUMED  TO  BE  NON  FLUCTUATING 
C**TARGET  RETURN  POWER  IS  UNITY 

COMMON/ ADPAR/DELT  *OAD.VMAX 

COMMON/PROPAR/NLOOK  »NSAMPP*NFFT  .NSCLK.NCINT 
COMMON/xCOOE/NSTAGE»lNIT(lO»*FEEDBK<lO> .CODE (1024) 
COMMON/SIGPAR/FOOP.NOELTA 
COMPLEX  X(l) 

DIMENSION  XR (4) 

INTEGER  CODE 
TPI*6. 283185307 
CALL  RAN0U<XR.4) 

PSI=XR(2)«TPI 

OMDOP=FD0P*TPI 

TIME=0.0 

NCYCLE=NLOOK/NSAMPP 
DO  10  IC=1 .NCYCLE 
DO  11  10s 1 .NSAMPP 
J= ( IC-1 ) *NSAMPP. ID 
ICOOE=ID-NDELTA 

IF < ICODE. LE*0>  ICODE=ICODE.NSAMPP 

ARG=TIME«OMOOP.PSI 

XI=FL0AT (CODE ( ICODE) ) *COS (ARG) 

XQ=FLOAT (COOE ( ICODE) ) *S1N ( ARG) 

X ( J) =X ( J) .CMPLX (XI »X0) 

11  TIMEsTImE.DELT 
10  CONTINUE 
RETURN 
END 
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SUBROUTINE  TAYWG <X*NLOOK> 

£«*#**••**#•*********«#**«<»•*«»*•»******«»«»*****»********«•#•***«**# 

C**IMPLEMENTS  TAYLOR  WEIGHTING  FOR  ANALOG  PROCESSOR 
C**X  IN  THE  COMPLEX  INPUT/OUTPUT  ARRAY 
C**NLOOK  IS  NUMBER  OF  SAMPLES  IN  X 

. . * . . . 

COMPLEX  X(l) 

OIMENSION  F (9) 

OATA  F/.462719*  ♦ 12681 6E-1 * •302744F-2*-* 178566E-2* •884107E-3* 
*-.38243?E-3«.121447E-3»-.4l7574E-5»-.249574E-4/ 

N=13500 

N2=N/2 

OEL ANG=6 • 283 1 8530 7/N 
DO  10  K=1*NL00K 
WEIGHT=1.0 
DO  20  M=1 » 9 
INCsM«(K-l-N2) 

ANGsDELANG«INC 

20  WEIGHT=WEIGHT*F (M> *COS (ANG) «2#0 
10  X(K>=X(k)«WEIGHT 
RETURN 
END 
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SUBROUTINE  WEIGHT <X»N> 

C»***************«********«*»**»************ 

c»*hamming  weighting  prior  to  fft  processing 
******♦**«»•«***»*«****•*»*• ••*****•**««* 
DIMENSION  X(l) 

COMPLEX  X 
N2=N/2 

PI=3.14J592654 

OELANG=Pl/FLOAT(N) 

ANGs-PI 
00  10  1=1 *N 
ANG=ANG*DELANG 
WIND=0.SA*0.46*C0S(ANG) 

X(n=X<I)*WIND 
10  CONTINUE 
RETURN 
END 
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